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ABSTRACT 

We compare the average star formation (SF) activity in X-ray selected AGN hosts with a mass-matched control sample of inactive** 
galaxies, including both star forming and quiescent sources, in the 0.5 < z < 2.5 redshift range. Recent observations carried out 
by PACS, the 60 - 210 pm photometric camera on board the Herschel Space Observatory, in GOODS-S, GOODS-N and COSMOS 
allow us to obtain an unbiased estimate of the far-IR luminosity, and hence of the SF properties, of the two samples. Accurate AGN 
host stellar mass estimates are obtained by decomposing their total emission into the stellar and the nuclear components. We report 
evidence of a higher average SF activity in AGN hosts with respect to the control sample of inactive galaxies. The level of SF 
enhancement is modest (~ 0.26 dex at ~ 3c confidence level) at low X-ray luminosities (L x < 10 43 ' 5 erg s~') and more pronounced 
(0.56 dex at > lOcr confidence level) in the hosts of luminous AGNs. However, when comparing to star forming galaxies only, AGN 
hosts are found broadly consistent with the locus of their 'main sequence'. We investigate the relative far-IR luminosity distributions 
of active and inactive galaxies, and find a higher fraction of PACS detected, hence normal and highly star forming systems among 
AGN hosts. Although different interpretations are possible, we explain our findings as a consequence of a twofold AGN growth path: 
faint AGNs evolve through secular processes, with instantaneous AGN accretion not tightly linked to the current total SF in the host 
galaxy, while the luminous AGNs co-evolve with their hosts through periods of enhanced AGN activity and star formation, possibly 
through major mergers. While an increased SF activity with respect to inactive galaxies of similar mass is expected in the latter, we 
interpret the modest SF offsets measured in low-Lx AGN hosts as either a) generated by non-synchronous accretion and SF histories 
in a merger scenario or b) due to possible connections between instantaneous SF and accretion that can be induced by smaller scale 
(non-major merger) mechanisms. Far-IR luminosity distributions favour the latter scenario. 
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1. Introduction 

Star formation (SF) and the phenomenon of Active Galactic 
Nuclei (AGNs) are two well-known and well studied aspects 
of galaxies. The former is related to the growth of stellar mass 
while the latter is tied to the growth of nuclear super-massive 
black holes (SMBHs). Both processes play primary roles in the 
formation and evolution of galaxies. 

Most star forming galaxies over most of the age of the 
Universe follow a strong correlation between their star forma- 
tion rate (SFR) and stellar mass, defining the so-called main 



Send offprint requests to: P. Santini, e-mail: 
paola. santini@oa-roma. inaf . it 

* Herschel is an ES A space observatory with science instruments pro- 
vided by European-led Principal Investigator consortia and with impor- 
tant participation from NASA. 

** Throughout the paper, the wording 'inactive'/'active' refers to 
galaxies lacking/showing nuclear activity (non-AGNs/AGNs), regard- 
less of their star formation rate. 
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sequence of star format ion. This tight relationship appears to 
hold up to z ~ 6 (e.g . iBrinchmann et aD l2004t [Noeske et al.l 
20071: lElbaz et al 1120071: fDaddi et al.ll2007bt ISantini et al.ll2009t 
Stark et al.1 120091) . While its normahzation is known to in- 
crease with redshift up to z ~ 2 - 3, an epoch which cor - 
responds to the peak of SF (e.g. IHopkins & Beacoml 120061) . 
and to flatten at higher redshifts, its slope is still debated 



iig_ 

(see discussion in IDunne et al.l 12009b IPa nnella et al. 2009; 



ISantini et alj|2009t iRodighiero et al J 1201 (It iKarim et all 1201 lb . 



Outliers from the main sequence are also observed, such as 
Ultra luminous Infrared Galaxies (ULIRGs, I Sanders & Mirabel 
19961) and Submillim e ter Galaxies (SMG s) (Chapma n et al 
2003b iPope etai1l2008b iTacconi et al.ll2008l) . These are charac 



terized by very large SFRs and specific SFRs (SFR per unit stel- 
lar mass). Studies based on CO emission interpret outliers as 
gas rich systems likely experiencing a different SF regime with 
respect to normal star forming galaxies on the main sequence: 
while the latter population is thought to form by smooth, secu- 
lar accretion with long duty cycles (e.g. iDekel et al.ll2009h . out- 
liers are believed to undergo short and intense starbursts occur- 
ring typically during merger^] or in dense n uclear star forming 
regions, which increase their SF^ efficiency dDaddi et aL 1201 Ot 
ITacconi et al1l2010tlGenzel et al.ll2010tlWuvts et alJl2oTll) . 



Observational evidence indicates that all massive galax- 
ies experience one or more phases of nuclear activity (e.g. 
iRichstone et ail 119981) during which their SMBHs grow. 
Several observations support a scenario of close connec- 
tion between AGNs and their hosts. Tight correlations are 
known to exist between SMBH mass and galaxy bulge 



cal Universe (e.g. Gebhardt et a 


. 20001 Ferrarese & Merrittl 


20001 Marconi & Huntl 12003b 


Kormendv & Benderl l2009t 


Giiltekin et al. 2009t iGraham et al. 


1201 lh and are also traced to 


higher redshift (e.g. iMerloni et al 


1 120101 iBennert et alJ 1201 lb. 



In addition, galaxies and SMBHs demonstrate similarities in 
their luminosity-dependent evolution with redshift: the space 
density of luminous AGNs peak at z ~ 2, while t hat of lower 
luminosity AGNs has it s maximum at z ~ 1 (e.g. Cowie et al 



2003 



2005 



Fioreetal. 2003 



a Franca et al.l 12005b iHasinger et al 



Brandt & Hasingeri l2005t iBongiorno et al] 120071 ). This 



so-called anti-hierarchical evolution is similar to the downsi zing 
behaviour of galaxy SF activity (see iFontanot et al.l 120091 and 
references therein). Finally, AGN growth seems to be mostly 
due to matter accretion from the host galaxies during the active 
phases of the AGN (although me rgers of black holes may be 
responsible for part of the growth. [Marconi et al. 2004). Bright 



ated with intense starburst events (e.g. 


lowan-Robinsonll 1 995t 


Barvainis & Ivisonl 20021 Omont et a 


J l2003t iPriddev et alJ 


20031 iPage et al.l 120041 IStevens et al 


12005b. Post-starburst 



QSOs. 



spectral signatures have been reported 
well as in lower luminosity l ocal AGNs (ICanalizo & Stockton 



lao et al 



well as in lower lumin osity local AGNs ICanalizo & S 
l20^ lKauffmann et al.ll2003l: iDavies et alll2007b . Inlsh 
(1201 Ob . we found a substantial level of star formation in both 
low and high luminosity AGNs, increasing with redshift at 
simi lar rate as the star formation in ina ctiv e galaxies (see 
also iLutz et all 120101: IMullanev et ai]l2012l and IMullanev et al] 
120 ioF Finally, a number of studies reported evidence for a 
correlation between AGN and far-infrared (FIR) luminosity, 
as well as between AGN luminosity and PAH feature strength 



1 Hereafter, unless otherwise stated, we refer to gas-rich major merg- 
ing when talking about 'mergers'. 



(iRowan-Robinsonl 119951 ISchweitzer et al.l 120061: iNetzer et all 
120071: ILutz et al.ll2008blBonfield et al.ll2o"TTl) . 

Although there is general agreement that relationships ex- 
ist between SMBHs and their host properties, the mechanism 
responsible for triggering the active phase is still debated. 
At high AGN luminosities, a process is at play which cor- 
relates AGN activity with star formation in the host galaxy. 
Major mergers are a g ood possibility according to early ob- 
servational studies (e . g. Stockton! "l982t iHeckmanet al.l 11984 
IStockton & Ridgwavl Il99lb iHutchings & Neffl Il992b and t o 
more recent theoretical predictions (e.g. ISpringel etal]|2005h : 
these transport large amounts of gas to the centre of the merging 
galaxies, feeding SMBHs and triggering intense SF episodes. 
Such a co-evolutio nary scenario is sup po rted by studies of 



s et al.l 
1200(1 



1988allbb ISarTders & Mirabel 



2001), which are all in ma- 



local ULIRGs (e.g. iSanders i 
11996b ICanalizo & Stockton! ! 
jor mergers, and whose exceptionally high SFR is often ac- 
companied by powerful AGN activity. Moreover, there is good 
evidence for a cl ose connection between AGN activit y and 
SF in SMGs (e.g. " 



Chapman et aT1l2005b I Alexander et al.ll2008b 



these studies may be subject to 
since such bright IR galaxies require, by con- 



ICoppin etai1l2008ir However, 
selection effects, 
struction, large gas amounts, which are supplied by major merg- 
ers. Alternative secular mechanisms of gas inflow, 'non-merger' 
mechanisms hereafter, have also been suggested as means to 
drive SF and trigger AGNs simultaneously. These mechanisms 
include minor mergers, disk instabilities and bars, supernova 
explosions or the inf all of recycl e d gas returned to the inter 



Genzel et al 
Ciotti et al.1 1201 Ob 



stellar medium (e g IWadal [2004b iMartinil 12004b Hogel 



2009; 



2008b IDekel et alJ 12009b iJohansson et al 

While major mergers are usually invoked 



2006 



uasars (e.g. Hopkins & Hernquistl l2009b 



to produce bright 

IVeilleux et al.l l2009bt iFiore et al.l 12012b . secular processes are 
asserted to be suffic ient to trigger the active nucleus in low lumi- 



nosity AGNs (e.g. IHopkins & Hernquistl 12009b IMullanev et al.1 



iity 

2012b. However, r ecent studies (e.g. ICisternas et al.l 1201 1 



Allevato et al.ll20llT) support evidence for major mergers not be- 



ing the leading triggering mechanism even in moderately lumi- 
nous AGNs (L x ~ 10 43 - 5 erg s" 1 ). 

AGNs are believed to play an important role in regulating 
their host galaxy's SF activity, and are therefore considered to 
be fundamental ingredients of theoretical models of galaxy for- 
mation and evolution. Nuclear emission is believed to be re- 
sponsible for gas heating and consequently SF suppression. This 
'negative AGN feedback' is necessary for theoretical models to 
reproduce many observations, such as t he strong suppression 



of SF in the most mas sive galaxies (e.g. i Di Matteo et a~l 
ISpringel et ail 12005b iBower et all 12006b ICroton et all 



2005 



2006 



iMenci et alj|2008t IHopkins & Elvisll2010b . Moreover, 'positive 
feedback' may also occu r, where AGN-driven winds induce SF 
in the host galaxy (e.g. iBegelman & Cioffil Il989t . ISiikl 12005b 
iFeain et al.ll2007T:ISilk & Normanll2009HElbaz et all[200% . 

For all the reasons described above, a detailed knowledge 
of the interplay between AGN accretion and star formation pro- 
cesses of the host galaxy is needed to reach a full understanding 
of galaxy formation and evolution. To shed light on their con- 
nection, we study the SF activity in X-ray selected AGN hosts 
from z ~ 0.5 to z ~ 2.5, i.e. in the range of redshifts over which 
most of the stella r and SMBH conte nt of the Universe has been 
put into place. In lShao et al. I J2010h we made the first effort in 
this direction, and we used Herschel SDP data in GOODS-N 
to study the dependence of host star formation on redshift and 
AGN luminosity. We found that the level of SF in AGN host 
galaxies increases with redshifts at a similar rate as in inactive 
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galaxies. Moreover, at each given redshift, we found no depen- 
dency between the host FIR luminosity (used as a proxy of the 
SFR) and the AGN luminosity up to Lx ~ 10 44 erg s 1 , while 
a corr elation is observed in brighter AGNs (Fig. 6 of Shao et al. 

We interpreted these results as reflecting the interplay of 
two paths of AGN-host co-evolution: very luminous AGNs are 
closely coupled to their host galaxy's growth by an evolutionary 
mechanism, for which one possible explanation is major merg- 
ing; on the other hand, no close coupling is observed between 
AGN accretion and total host SFR in low luminosity AGNs, re- 
flecting a secular evolution. The level of this secular star for- 
mation increasingly dominates over the correlation at increasing 
redshift. 

This study expands upon our previous work. We take the 
next step and systematically compare the SF properties of AGN 
hosts to inactive galaxies of similar stellar mass in order to ac- 
count for the important covariances of stellar mass with SF. For 
this purpose, accurate estimates of the host stellar masses were 
obtained by means of a specific decomposition technique de- 
veloped to separate stellar and nuclear emission. Measuring the 
SFR in AGN hosts is not trivial, since nuclear emission can 
frequently outshine both the ultra-violet (UV) and mid-infrared 
(MIR) emission from young stars. However, the FIR contin- 
uum is shown to be dominated by the host galaxy emission for 
all but the most extreme AGN-dominated sy stems dNetzer et alJ 
120071: iLutz e"tai1 l2008l: iMullanev et al1l201 ll) and to be a proxy 
of its SF activity ISchweitzer et al.ll2006l: ILutz et al.ll2008l) (see 
also S ect. 13.21 ). The Herschel Space Observatory ( Pilbratt et alJ 
12011 allows us to probe the FIR nature of galaxies with greater 
depth and angular resolution than possible with previous FIR 
and sub-mm facilities. The Photodetector A rray Camera and 
Spectrometer (PACS, iPoglitsch etail 1201 Oh . with its 60-210 
fim window, is able to detect the FIR emission of dust heated by 
UV photons from young stars. Herschel allows a clean measure- 
ment of the SF activity in a wavelength range that is relatively 
unaffected by AGN contamination or uncertain UV dust correc- 
tions. In addition, PACS studies are not restricted to the modest 
sample sizes and severe flux-limited selection of MIR spectro- 
scopic campaigns. 

The structure of the paper is the following: after introducing 
the sample and the selections applied to it in Sect|2] as well as the 
methodology used to infer stellar masses and FIR luminosities 
(and hence SFRs) in Sect.[3j we present our results in Sect. [4] In 
Sect. [5] we discuss the implications of this analysis and suggest 
possible interpretations. Finally, Sect. [6] summarizes the major 
steps and results of this work. 

We adopt a H = 70 km s _1 Mpc~', Q. M = 0.3 and Q A = 0.7 
cosmology and assume a Salpeter IMF. 

2. Data set 

This work exploits the excellent multi-wavelength coverage 
available in three deep extragalactic fields: GOODS-South, 
GOODS-North and COSMOS. The depth of the two GOODS 
fields in the X-ray, FIR and several other bands is essential 
to probe faint and high redshift galaxies and AGNs. The shal- 
lower, wide-area COSMOS field provides us with good statistics 
among massive and bright sources, which are rare in the GOODS 
fields. 



2.1. Far-infrared (FIR) data 

The far-infrared d ata used in this work were collected by the 
PACS instrument dPoglitsch et al.ll20loh on board the Herschel 



Space Observatory dPilbratt et al.l l2010h. as part of the PACS 
Evolutionary Probe (PEP. lLutz et al.l 201 ll) survey. Observations 
of GOODS-S were carried out in all three PACS bands (70, 100 
and 160 fim), while the other two fields were only observed 
in the two long wavelength bands. We made use of PACS cat- 
alogues extracted based on the prior knowledge of the posi- 
tions and fluxes of source s detected in deep archival MIPS 24 
m imaging in these fields dMagnelli et al.l2009tlLe Floc'h et all 
20091) . This allows us to accurately deblend PACS sources in im- 
ages characterized by a large PSF, especially in crowded fields, 
and permits us to improve the completeness and reduce the num- 
ber of spurious sources at faint levels compared to a blind extrac- 
tion. Fluxes in GOODS-S reach 1.1, 1.2 and 2.4 mJy at 3cr in the 
70, 100 and 160 fim bands respectively. In the other two fields 
the 3cr flux limit at 100/160 /mi is 3.0/5.7 mJy in GOODS-N 
and 5.0/10.2 mJy in COSMOS. We consider sources below these 
limits as undetected by PACS. Detailed information on the PEP 
survey, observe d fields, da t a proc essing and source extraction 
may be found in lLutz et al.l (1201 ll) . 

The technique used to associate PACS fluxes with optical 
counterparts (i.e. PACS-to-24 fim and 24 /im-to-IRAC and opti- 
cal association based on IRAC positions) could in principle in- 
troduce some biases when applied to the specific topic of this 
work, namely the comparison of FIR properties of AGN hosts 
and inactive galaxi es. The first reason i s that, given their strong 
MIR emission re. g. iDaddi eTaIll2007aL [Fiore et al.ll2008l) . AGN 
hosts are more represented in the 24 fim prior list compared to 
inactive galaxies (for a given FIR host's luminosity, they have a 
larger probability of being detected at 24 fim). We might there- 
fore associate FIR flux to an AGN host rather than an inactive 
galaxy if the latter is not represented in the prior list. However, 
the number of PACS sources wi th no 24 fim counterpart is very 
small (of the order of < 2%, iMagdis et all 12011b ILutz et al.l 
l20TTl) . The second source of possible bias is that, in case of se- 
vere blending in the PACS images, the flux is associated with 
the brighter MIR counterpart, which, again, is more likely to 
be an AGN (in case of one inactive and one active galaxy of 
similar host's luminosity falling in the same beam). However, 
such blending issues are quite rare (only ~ 4% of PACS sources 
would be affected in the reddest PACS band in GOODS while 
the fraction would be negligible in COSMOS). We performed 
two tests to understand the degree of these possible biases on 
our results. First, we repeated our an alysis only using PACS cat- 
alogues obtained by blind extraction (ILutz et al.l201 ll) . Then, we 
performed a prior extraction without using a MIR flux weighting 
scheme and re-ran our analysis. In both cases we found no sig- 
nificant difference in our results, given the very similar FIR flux 
distributions. These tests guarantee that the results presented in 
this work are not caused by artefacts occurring during source 
extraction and/or optical associations. 

2.2. Ancillary data 

In order to infer redshifts and stellar masses needed for this 
study, we complement PACS observations with multiwavelength 
photometric catalogues. 

In the two GOODS fields we used catalogues which provide 
photometry in 14 bands, from the U -band to 8 fim. For GOODS- 
S we used the updated GO ODS-MUSIC catalogue dSantini et al.l 
l2009HGrazian et al.l2006l) . Since the latter covers an area smaller 
than the entire GOODS-S field, for AGNs lying outside the 
GOODS-MUSIC footprint we us ed the photometric informa- 
tion compiled by iLuo et all d2010l) from other public photomet- 
ric datasets in the GOODS-S field. We verified that systematic 
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differences in the photometric catalogues were not generating 
serious inconsistencies in the stellar mass estimates (Sect. 13. U 
by comparing masse s of sources in common to GOODS-MUSIC 
and lLuo et al]d2010l) derived from both sets of multi-wavelength 
photometry. For GOODS-N w e used the catalogue com piled by 
the PEP Team and described in lBerta et all d2010l [20TTb . 

For the COSMOS field we used the 
multiwavelength public catalog available at 
http://irsa.ipac.caltech.edu/data/C OSMOS/tables/pho tometry/ 
The data reduction is described in Capa ketal](l2007t) . However 
the new catalog uses better algorithms for source detection and 
photometry measurements. Still, by comparing the COSMOS 
colour-mass diagram to that of GOODS, we found a sus- 
piciously large overabundance of faint red galaxies in the 
COSMOS catalogue, especially at high redshifts. The Subaru 
/-band photometry, which is the selection basis for this cata- 
logue, shows an extended tail towards fainter values for faint 
objects when compared to the ACS /-band. Extremely low 
Subaru / fluxes can then give rise to fake breaks in the observed 
spectral shape which may be confused with red and massive 
galaxy features. The requirement of at least a detection in the 
ACS /-band solved this problem, so we restricted both the 
AGN and the non-AGN sample to the area covered by ACS in 
order to ensure good quality photometry for all object in the 
COSMOS field. The COSM OS optical catalog is suppl emented 
with K photom etry from iMcCracken et all (1201 Ol) IRAC 
photometry from ISanders et alJ d2007l) and lllbert etalf d2009t) . 
and 24 //m photometry from iLe Floc'h et al. ( 20091) . As far 
as the control sample of inactive galaxies is concerned (see 
Sect. 12.41 1. instead of using the entire (more than two million 
sources) dataset, we randomly extracted from the central and 
most covered area a catalogue listing ~ 65000 ACS /-selected 
sources. 

All the catalogues are supplemented with either spectro- 
scopic redshifts, or photometric redshifts when spectroscopic 
ones are unavailable. Photometric redshifts in GOODS-S were 
computed by fitting the multiwavelength photometry to the 
PEG ASE 2.0 templates (iFioc & Rocca-Volmerangdll997l see 
details in Graz ian et al 1 120061 update d as inlSantini et al.ll2009h 



2CX)lj 



and by adopting the EAZY code ( iBrammer et al. 
GOODS-N and COSMOS (details in iBertaet all 1201 lh . The 
fraction of total outliers, defined as objects having |A(z)|/(l + 
Zspec) > 0.2, is 5.3%, 5.8% and 0.7% in GOODS-S, GOODS- 
N and COSMOS, respectively. Once these outliers have been 
removed, the average |A(z)|/(l + z S pec) is equal to 0.04 + 0.04, 
0.03 +0.03 and 0.01 +0.02, respectively in the three fields. We re- 
fer to the papers cited above for more detailed information about 
the photometric redshift estimate and their quality. 



2.3. X-ray information 

Detailed X-ray point source catalogues were compiled for all 
three fields. These provide rest-frame X-ray luminosities and ob- 
scuring column densities Nh derived from X-ray spectral anal- 
ysis, source classification, optical associations as well as, where 
needed, photometric redshifts specifically suited for AGN hosts, 
based on optimised galaxy, AGN and galaxy/ AGN hybrid tem- 
plates. 

For GOODS-S we used the lLuo etaf](l2008h catalogue from 
the 2 Msec Chandra Deep Field-South imaging program. Cross- 
matches, photojnetric_redshifts and X-ray properties are pro- 
vided bv lLuo et ai]d2010l) . 

For GOODS-N we used the 2 Msec C handra Deep Field- 
North catalogue of lAlexander et alJ d2003h and an updated ver- 



sion of the classificat ion into AGNs and other type s of X-ray 
sourc es compiled b y iBauer et al] (120041) (see also IShao et all 
l2010h . We refer to IBauer et all (2004) for details on spectro- 
scopic and photometric redshifts and X-ray properties. 

Finally, for CO SMOS, we used the XM M-COSMOS cata- 
logue compiled by ICappelluti et al] d2009l) and the X-ray opti- 
cal ass ociations and derived properties presented in lBrusa et al] 
(l2010h in its updated version, which includes new redshifts 
from ongoing (mostly DEIMOS/KECK) spectroscopic cam- 
paigns and few changes in the published redshifts and/or spec- 
troscopic classifications from a re-analysis of the data. X-ray 
luminosities and obscuring column densities are derived by 
Mainieri et al. 1(120^12011 . Photometric redshifts are provided 
bv lSalvato et al ■I d2009h~ As mentioned in Sect. 12.21 we restricted 
ourselves to X-ray selected AGNs that lie in the area covered by 
ACS /-band, which are roughly 82% of the full XMM-COSMOS 
catalogue. 

A sense of the differing depths of the X-ray data in the three 
fields may be gained from Fig.Q] which shows the distribution of 
absorption-corrected rest-frame 2-10 keV luminosities {Lx here- 
after) as a function of redshift for the final AGN sample used 
in this analysis (see Sect. 12.41 for the selections applied to the 
total sample). GOODS (South+North) and COSMOS can be re- 
garded as complementary in terms of their X-ray properties: the 
deep GOODS fields allow a good sampling of faint AGNs, while 
in COSMOS we probe the rarest and most powerful sources. 

2.4. Sample selection 



As discussed in Sect. 13.11 we use the classifications of X-ray 
sources into optical Type 1 and Type 2 as a prior in the stel- 
lar mass estimates of AGN hosts. As a first step, we selected 
AGN host galaxies based on the classification provided by the 
X-ray catalogues. This classification is limited to optical spec- 
troscopic criteria for most of the sources in the GOODS fields 
and for the entire COSMOS sample. In order to increase the 
statistics and to avoid spectroscopic biases, we expanded the 
AGN sample by including all X-ray sources with Lx > 10 42 
erg s _1 even if not classified as AGN. We considered them as 
Type 2 in GOODS, being highly unlikely that Type 1 AGNs 
were missed by the dedicated spectroscopic campaigns that tar- 
geted AGNs in these deep fields. In COSMOS we instead used 
the clas sification pro v ided b y the photometric redshift proce- 
dure of ISalvato et al.l d2009h, which has been demonstrated to 
be hi ghly trustw orthy (e.g. Cappelluti et al.l 120091: iLusso et all 
12010b . Following ITusso et al.l d201oir we considered Type 1 all 
the objects best-fitted by a QSO or hybrid Type 1 - galaxy tem- 
plate (model SED > 18 in lSalvato et~aI1l2009l) . and Type 2 all 
the remaining objects. 

In order to derive reliable estimates of the FIR luminosity 
(Sect. l3~!2l . we only used sources in regions of our PACS images 
where the 100 fim coverage (i.e. integration time) is at least half 
of its value at the centre of the image. In the redshift range 0.5 - 
2.5 studied in this work, this coverage restriction removes ~ 2 1 % 
of AGNs in GOODS-S and ~ 12% in GOODS-N. We divided 
our sample into three mass bins (log M[MO] = 9 - 10, 10 - 11 
and 11-12, see Sect. l3.1l for the stellar mass estimate) and three 
redshift intervals (z = 0.5-0.8,0.8-1.5 and 1.5 -2.5) which are 
narrow enough to keep the K-correction in the FIR at reasonable 
values within any interval. As explained in Sect. 13.21 up to z ~ 
2.5 PACS is able to probe a rest-frame window which traces the 
SFR in both active and inactive galaxies. 

A key aim of this work is to compare AGN hosts with in- 
active galaxies of similar mass. For that purpose, we selected 



4 



P. Santini et al.: Enhanced SFR in AGN hosts from PEP-Herschel observations 




0.5 1 1.5 2 2.5 
Redshift 



"i — 1 — i — 1 — i — 1 — i — 1 — r 



_L 



GOODS-N J 

J i I i L 



0.5 1 1.5 2 2.5 
Redshift 



"i — 1 — i — 1 — i — 1 — i — 1 — r 



. I. .-. .. ■ 



COSMOS _ 

J i I i L 



0.5 1 1.5 2 2.5 
Redshift 



Fig. 1. Intrinsic X-ray luminosities at 2-10 keV {Lx) as a function of redshift in the three fields for the final AGN sample used in the 
analysis (Sect. l2~4t . Horizontal red lines show the thresholds used to split the sample into high- and low-L^ subsamples (Sect. [570 . 



a reference sample of galaxies in the same PACS fields and in 
the same redshift and mass intervals as the AGN hosts' sam- 
ple. We excluded X-ray selected AGNs and Galactic stars, as 
well as all sources with 2-10 keV X-ray luminosity larger than 
10 42 erg s _1 , even if they were not explicitly classified as AGNs. 
Although this degree of X-ray emission mostly traces AGN ac- 
tivity, it may also ar ise from intense star formation at the level 
of a few 10 2 M o /vr jRanaUi et al] 120031: lAlexander et alJl2005t 
lLaird et alJ20"05ll201oir Removal of these sources may therefore 
bias the control sample towards a lower average SFR by remov- 
ing the most strongly SF systems. Such a bias would at the same 
time push the average SFR high for the AGN sample, if these 
intense starbursts were mistaken for AGNs. We tested the im- 
portance of strong starbursts among the low-L^ X-ray sources in 
two ways. Since such strongly star forming galaxies are certain 
to be detected with PACS and lie in the exponential tail of the 
population of SF galaxies at any given redshift, we first searched 
for strong variations in the PACS detection rates with Lx among 
low luminosity X-ray sources in GOODS-S. None were found. 
We also repeated the entire analysis of this paper after including 
X-ray sources with no AGN classification in the control sample 
instead. Our results were unchanged. We conclude that including 
or excluding these very strong starbursts has a negligible effect 
on our reference catalogue and does not affect our results to any 
degree. 

A sizeable fraction of very MIR bright galaxies, usually se- 
lected from their 24 fim excess, are likely to be highly obscured 
AGNs, whic h are too dust enshrouded to be detected by X-ray 
surveys (e.g. Daddi et al.l feOOTat iGeorgantopoulos et alj l2008t 
iFiore et al.ll2008l 120091) . Removing these objects as AGN can- 
didates would bias our FIR luminosity estimate, sin ce some of 
these g alaxies are also actively forming s tars (e.g. iLutz et al 



tnese gaiaxies are aiso activeiy torming stars (e.g. Lutz et ; 
2005t ISaiina et alJ l2008t iDev et al.ll200a iDonlev et al.l 1201 



s 



Georgakakis et al.|l2010h . However, no more than 0.2% of the 



sources in our inactiv e samp le satisfy the selection criterion sug- 
gested bv lFiore et alj (120081) . We kept these sources in our con- 
trol sample, and we checked that the results presented below re- 
main unchanged when they are removed. 

For consistency with the COSMOS sample, which is re- 
stricted to /-detected sources in order to ensure good photomet- 
ric quality (Sect. l2~!2l . we applied the same optical selection to all 



the samples considered, namely all AGNs and all non-AGNs in 
all the fields. The /-band filters used in GOODS and COSMOS 
are slightly different, with the latter having 20% more power on 
the red side. However, given the good spectral sampling in the 
optical range, we do not expect significant bias effects from this 
small difference. 

The /-band selection may in principle introduce a bias in 
the analysis, because, for z > 1 galaxies, the observed /-band 
moves to UV rest-frame wavelengths. Dust obscured star form- 
ing and quiescent galaxies might therefore not enter the /-band 
cut and be missed in the analysis, and it is not clear whether and 
at what level this affects our comparison. However, we checked 
that, given the depth of the GOODS fields and the correlation 
between the /-band emission and the stellar mass, the samples 
extracted from these fields are insignificantly affected by the 
/-band cut in the mass range studied. The same does not hold 
for the shallower COSMOS field. However, the similarities ob- 
served between GOODS fields and COSMOS (see Sect. [57Tb 
make us confident that our analysis is unbiased by this effect. 

As we show in the next section, AGNs are typi- 



cally hosted by high mass galaxies (e.g. Bundv et al. _ 2008 ; 



Alonso-Herrero et al.l l2008t iBrusa et"aU I2009T 



ISilverman et al.l 

l2009bt IXue et al.l l2010t iMainieri et al.1 1201 ll) . The mass bins 
used in this analysis are large enough that a difference in the 
mass distributions of AGNs and control galaxies, even within a 
given bin, could add some bias to our estimates of mean FIR 
properties. We therefore restricted ourselves to a mass-matched 
reference sample of control galaxies, to take out any possible 
covariance between SFR and stellar mass due to differing mass 
distributions between active and inactive galaxies. We matched 
each AGN with 6 non-AGNs with a mass within +0.2 dex of 
the AGN host mass. Since AGNs typically reside in massive 
galaxies, we might run out of comparison galaxies at the high 
mass end, especially at high redshift. Under this circumstance, 
we allowed control galaxies to be picked more than once and as- 
sociated to these sources the appropriate weight. This way, we 
only lose one AGN in GOODS-S (due to the lack of inactive 
galaxies of similar mass). The mass tolerance and the number of 
matches were chosen to maximize the final control sample and 
at the same time not have too many repeated galaxies among the 
control sample (4% in the worst case, typically < 2%). 
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We also checked whether the mass-matched samples are af- 
fected by any redshift effect within each redshift bin. We con- 
sidered the difference between the average values of the red- 
shift distributions of AGNs and mass-matched non-AGN in the 
10 9 - 1O 12 M mass interval and in each redshift bin. The differ- 
ence A(z) = (z) n0 nAGN ~ (z)agn is of the order of ~ 1(T 3 - 1(T 2 
and consistent with zero (errors are computed by assuming nor- 
mal distributions) in most of the bins. The only redshift intervals 
where A(z) is significant are the high-z bin in COSMOS and the 
low-z and high-z bins in GOODS-N. The latter bin shows the 
largest disagreement between non-AGN and AGN redshift dis- 
tributions (0.18 + 0.04). However, as we will show in Sect. |4j 
this mild redshift inconsistency, with AGNs lying at an average 
redshift lower than non-AGNs, even improves the robustness of 
our results. 

To summarize, both AGN and mass-matched non-AGN sam- 
ples are characterized by: 

- /-band detection; 

- good PACS coverage; 

- 0.5<z<2.5; 

- 10 9 < M[M Q ] < 10 12 . 

Within these requirements, the two samples are selected accord- 
ing to the following criteria: 

• AGNs: 

- detected in X-ray; 

- classified as AGNs from optical spectroscopy; 

- additionally includes unclassified X-ray sources if Lx > 
10 42 erg s , considering them as Type 2 in GOODS and 
using the classification from the photometric redshift fit- 
ting with hybrid templates in COSMOS; 

• non-AGNs: 

- excluded all AGN-classified sources; 

- excluded Galactic stars. 

After applying all the selections above, we ended up with 
136, 159 and 1052 AGNs in GOODS-S, GOODS-N and 
COSMOS respectively in the stellar mass range of 10 9 - 10 12 M Q . 
Of these, 19, 11 and 439 are Type 1 and 117, 148 and 613 are 
Type 2. The non-AGN reference samples are, by definition, ~ 6 
times larger than the active ones. The fraction of spectroscopic 
AGNs is 74%, 31% and 57%, and the fraction of spectroscopic 
control inactive galaxies is 47%, 43% and 5%, respectively in 
the three fields. By performing, as a control check, the analy- 
sis described below on purely spectroscopic sources, we verified 
that the difference in the spectroscopic fractions among the var- 
ious fields and between AGN and non-AGN samples does not 
bias our results. 

The statistics in the different redshift and mass bins for all 
three fields is summarized in Tab.Q] 

3. Method 

In this section we describe how we derived the fundamental in- 
gredients of this analysis, namely stellar masses and far-infrared 
luminosities. Results presented in Sect. 13.11 apply to the entire 
sample of AGNs and non-AGNs in the redshift range 0.5 - 2.5. 
The selection based on PACS coverage (Sect. l2~4l i was neglected 
since it does not affect stellar mass estimates. 



3.1. Stellar masses 

Stellar masses are most robustly estimated by fitting observed 
photometry with a library of stellar synthetic templates (e.g. 



B 19 




Fig. 2. Observed and best-fit SEDs of a Sy 1 galaxy (bottom) and 
a Compton Thin Sy2 (top), both lying in the GOODS-S field. 
Left panels show the fit with pure stellar templates, while the 
right ones show the decomposition between the stellar (blue) and 
nuclear (red) components. In the latter, the black lines represent 
the total emission (stars + AGN). The inferred stellar masses are 
reported in each panel. 



iFontana et alJ 1200*61) . For our set of inactive gal axies, masses 



were derived using the same procedure as in ISan tini et al 
( 2009 ). performing a x 1 minimization of iBruzual & Chariot 
(l2003h synthetic models, assuming a Salpeter IMF and param- 
eterizing the star formation histories (SFH) as exponentially de- 
clining laws. Each band is weighted with the inverse of the pho- 
tometric uncertainty. l<x errors determined by photometric un- 
certainties were computed b y accounting for all the s olutions 
within xleduced mm + Since IBruzual & Charlol (l2003h models 
do not include emission from dust reprocessing, we fitted the 
observed fluxes up to 5.5 fim rest-frame, fixing the redshift to 
the photometric or spectroscopic one, where available. 

Of course, this procedure is only reliable under the assump- 
tion that the bulk of the light emitted by the galaxy comes from 
its stellar component. This assumption is no longer trustworthy 
in the case of AGN host galaxies. In Type 2 AGNs, near- and 
mid-infrared bands are boosted by the dust torus emission, typi- 
cally biasing the mass estimate towards higher values (although 
at some level star formation can also affect the 3.6 //m band, 
iJosephet all 1 1984 iGlass & Moorwood [19851) . The situation is 
even more critical for Type 1 AGNs, whose Spectral Energy 
Distribution (SED) is often completely dominated by the nuclear 
emission at even shorter wavelengths. To account for the nuclear 
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Table 1. Statistics of AGNs and inactive sample. 



GOODS-South + GOODS-North 








AGN hosts 






inactive galaxies 








10" - 10 10 M Q 


10 10 - 10"M Q 


10" - 10 12 M o 


10" - 10 10 M G 


10 10 - 1O U M 


10" - 1O 12 M 


z = 0.5 
z = 0.8 
z= 1.5 


-0.8 
- 1.5 
-2.5 


9 (33.3%) 
19(15.8%) 
5 (0.0%) 


34 (20.6%) 
70 (18.6%) 
59 (6.8%) 


14 (35.7%) 
49 (30.6%) 
36 (33.3%) 


49 (4.1%) 
114 (0.0%) 
34 (0.0%) 


228 (18.9%) 
458 (11.6%) 
377 (4.2%) 


59 (15.3%) 
230 (12.6%) 
168 (17.3%) 


COSMOS 








AGN hosts 






inactive galaxies 








10' - 10 I0 M G 


10 10 - 10"M o 


10" - 10 12 M o 


10 9 - 10 10 M Q 


10 10 - 10"M o 


10" - 10 I2 M G 


z = 0.5 
z = 0.8 
z= 1.5 


-0.8 
- 1.5 
-2.5 


19 (0.0%) 
15 (6.7%) 
3 (33.3%) 


75 (8.0%) 
249 (4.0%) 
170 (4.1%) 


89 (14.6%) 
245 (13.9%) 
187 (9.1%) 


114 (0.0%) 
117(0.0%) 
24 (0.0%) 


474 (4.0%) 
1689 (1.2%) 
1126 (0.3%) 


507 (5.7%) 
1071 (4.4%) 
592 (1.9%) 



Notes. Number of AGNs and inactive galaxies in each redshift (rows) and mass (columns) bin in the different fields. Numbers in brackets show 
the fraction of sources fully detected by PACS (i.e. detected in the two bands used to infer vL v (60yum) as explained in Sect. l3.2t . 



emission and to achieve a reliable estimate of the stellar mass 
in AGN host galaxies, we created a special template lib r ary by 
combining the stellar emission from iBruzual & Chariot! d2003h 
model with the nuclear emission described by the Silva et all 
d2004l) templates. Once the total observed SEE) is decomposed 
into a stellar and a nuclear component, the stellar mass, as well 
as the other stellar parameters, are inferred from the stellar com- 
ponent only. 

The ISilva et al.l (120041) templates consist of four (one un- 
obscured, three affected by different levels of obscuration) nu- 
clear SEDs of Seyfert galaxies, for which the stellar contribu- 
tion was removed, normalized to the hard X-ray (2-10 keV) in- 
trinsic luminosity and averaged within bins of absorbing Nh- 
We averaged the SEDs with 10 22 < N H < 10 23 cirT 2 and 
10 23 < N H < 10 24 cirT 2 to yield one Syl and two Sy2 
(Sy2 C 77im and Sy2 C Thick with respectively Nh < and > 10 24 
cm ' 2 ) templates. Before com bining each of these SEDs with 
the lBruzual & Charlotl(l20 03) library, to mi mic real galaxy ernis- 
sion we reddened the Syl template using a Calz etti et al.1 (2000) 
law with E(B-V) va lues between and 0.3 dHopkins et al.ll2004t 
iMerloni et al.l2010l) and step of 0. 1 . With the aim of reducing the 
degeneracies given by the large number of free parameters (stel- 
lar parameters, type of nuclear template, normalization of the 
nuclear contribution vs the stellar one, reddening), we created 
three different libraries (one for each AGN template) and used 
prior knowledge of X-ray properties (classification and Nh) to 
fit each AGN host galaxy with the appropriate two-components 
library: objects classified as Type 1 were fitted with the stel- 
lar+Syl library, Type 2 AGNs with N H < 10 24 «7r 2 were fit- 
ted with the stellar+Sy2c77 lln library while for heavily obscured 
AGNs (N H > 10 24 cw -2 ) we adopted the stellar+Sy2 C7 -/„ CJ t li- 
brary. 

We used redshifts from the X-ray catalogues, which provide 
updated/new spectroscopy as well as photometric redshifts ob- 
tained using specifically suited AGN templates (see references 
in Sect.l2~3l. 



In order to make the fitting procedure computationally fea- 
sible, we resorted to a library for the stellar component with a 
reduced number of p arameters. We kept t he same range of stel- 
lar parameters as i n I Santini et alj d2009l) (see also Table 1 in 
iFontana et al.l2004T) . but adopted a coarser grid, and we assumed 



a lCalzetti et alj d2000h extinction law. We checked that the use of 
this coarse grid did not introduce significant differences with re- 
spect to the finer grid used in previous works on the GOODS-S 
field. The distribution of (M coarse grid - M fine gri d)/M coarse grid for 
GOODS-S non-AGNs in the redshift range used in this work is 
sharply peaked around 0, with a very low median/mean value of 
~ 7 ■ 10~ 3 /1 .2 • 10~ 2 and a semi-interquartile range/standard devi- 
ation of ~ 0.09/0.28; the ratio of the two stellar mass estimates 
is higher than a factor of 2 in 3.6% of objects. For consistency 
with the AGN sample we apply the same coarse parameter grid 
to mass estimates for the non-AGN reference sample as well. 

We are aware that the exponentially declining models (i.e. t- 
models) adopted for building the stellar c omponent are likely a n 
oversimpl ified descr i ption f or SFH (e.g.. lMaraston et ai1l2010l) . 
However, iLee et al] d2010l) showed that the resulting stellar 
masses can still be considered robust because of a combina- 
tion of effects in the estimate of the galaxy star formation rates 
and ages. Moreover, they are widely used in the literature and 
allow an easy comparison with previous works. Nonetheless, 
their application to active galaxies can be even more problem- 
atic, since we expect that, at least in a fraction of them, a burst 
of star formation is triggered at some point during their evolu- 
tion, making the declining tail of the SFH an even less accu- 
rate description. We therefore repeated the mass estimate by us- 
ing exponentially increas ing star formation his tories (i.e. the so- 
called inverted T-models. lMaraston et al.ll2~010t) . which represent 
the opposite approximation and may resemble better the SFH in 
AGN hosts galaxies. We found very low differences in the best- 
fit stellar masses in GOODS-S: the median/mean value of (M T - 
Mi nverte d T )/M T is equal to 0.07/0.03 with a semi-interquartile 
range of 0.14 and a standard deviation of 0.31 (these numbers 
refer to the 0.5 - 2.5 redshift range). The effect of the SFH is 
slightly larger in COSMOS, where many more bright AGNs are 
detected: the median/mean value of {M T -M mverte d T )/M T is equal 
to 0.11/0.11 with a semi-interquartile range of 0.18 and a stan- 
dard deviation of 0.35 (again for the 0.5 - 2.5 redshift range). 
However, we can consider the star formation history parame- 
terization to have an effect on the stellar mass estimate which 
is lower than its intrinsic lcr uncertainty in both fields (see be- 
low). Moreover, T-models are preferred over inverted r-models 
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Fig. 3. Normalized distributions of the relative difference be- 
tween the stellar mass estimates when using pure stellar tem- 
plates (M) and the decomposition method (M dec ). Blue and red 
histograms represent the total sample of Type 1 and Type 2, 
respectively. We report the mean/median values and standard 
deviation/semi-interquartile range of each sample. 



from an analysis of the x 2 > in agreement with the results of 
iRosario et all (1201 ll) . 

The adoption of the widely used iBruzual & Chariot! (|2003) 
library instead of i ts more update d version including TP-A GB 
stars contribution dBruzua ll 120071 see also lMarastonll2005h al- 
lows an easier comparison with previous works. The improved 
TP-AGB stars treatment predicts a larger rest-frame A'-band flux, 
and hence slightly lower stellar mass estimates. By comparing 
the stellar masses obtained with the 2003 and the 2007 libraries, 
we found that the latter underpredicts stellar masses by 0. 12 dex 
on average, although the scatter between the two estimates is 
as large as 0.17 d ex and there is no c lear trend with redshift or 
stellar mass itself (ISantini et al.ll2.Q12h . However, AGNs and in- 
active galaxies should be affected by the same possible system- 
atics, and since we aim at performing an unbiased comparison 
between the two populations, the analysis is internally consis- 
tent. 

Figure [2] shows an example of the SED decomposition 
method described above. It presents a Sy 1 galaxy (bottom) and a 
Compton Thin Sy2 (top). Left panels show the fit with pure stel- 
lar templates, while the right ones show the decomposition be- 
tween the stellar (blue) and nuclear (red) components (the black 
curves showing the total best-fit SED). From these plots it is ev- 
ident that in these objects pure stellar templates can not properly 
reproduce the emission at long wavelength (and even at short 
wavelengths in the most powerful AGN hosts). 

The typical Type 1 SED is usually highly dominated by the 
nuclear component, and sometimes almost completely outshone 
by it. Disentangling the stellar and the nuclear components in 
Type 1 AGNs is therefore very complex, and the derived stel- 
lar masses are affected by large uncertainties. The average lcr 
relative uncertainty (+ standard deviation) in the stellar mass 
estimate AM/M, where AM - (M max - M min )/2 for inactive 
galaxies and for Type 2 is 0.16 + 0.10 and 0.18 + 0.11 respec- 
tively in GOODS-S, 0.19 + 0.10 and 0.22 + 0.09 respectively 
in GOODS-N, and 0.30 + 0.13 and 0.30 ± 0.11 respectively 
in COSMOS. The average lcr relative uncertainty in the stel- 
lar mass estimate for Type 1 is significantly larger, being equal 



to 0.49 + 0.22, 0.52 + 0.40, and 0.64 + 0.35, respectively, in 
GOODS-S, GOODS-N, and COSMOS. (These numbers refer 
to AGNs and mass-matched non-AGNs in the 0.5 - 2.5 redshift 
range). Given these caveats, we checked our results against pos- 
sible biases introduced by uncertain estimates of the stellar mass 
in Type 1 AGNs. By repeating the analysis described in Section 
4, but only on Type 2 AGNs, we verified that the principal con- 
clusions of this work are unaffected by these large uncertainties. 
We a lso compared our res ults with the sample of Type 1 AGNs 
from lMerloni et all (120101) . and we found good agreement when 
the same x 1 minimization method is used, despite the different 
AGN template used (see also Bongiorno et al. in prep.). 

In order to characterize the improvements in the stellar 
mass estimates achieved through the decomposition technique, 
we plotted in Fig. [3] the distribution of the relative difference 
(M - M dec )/M, where M is the stellar mass computed from pure 
stellar templates and M dec is the output of the decomposition. 
The red histogram refers to Type 2 AGN hosts, the blue one 
to Type 1. When only Type 2 hosts are considered, the mean 
and median values for (M - M dec )/M are consistent with zero, 
and only in 1.3% of objects the difference in the stellar masses 
is larger than a factor of 2. However, the shift becomes more 
severe for Type 1 AGNs, modulo the uncertainties discussed 
above. Although the distribution goes in the direction of recov- 
ering lower stellar masses through decomposition (the median 
assumes a positive value, although still consistent with zero), the 
opposite is also found: in some cases the pure stellar fits might 
favour bluer, younger and lower mass objects in the attempt of 
reproducing the total (stellar + nuclear) light. The distribution of 
(M - M dec )/M in Type 1 AGN hosts is broad, and the difference 
in the stellar masses is larger than a factor of 2 in ~ 29% of the 
objects: the decomposition procedure seems to be essential to 
recover a reliable stellar mass in these sources. 

We also explored different ways to estimate the stellar mass 
in AGN hosts. The first one consists in fitting observed fluxes 
with pure stellar templates only up to 2 yum rest-frame (M2 ltm ), 
ignoring the near- and mid-infrared bands typically boosted by 
the AGN emission. As a second method (Magn subtraction), we 
used the X-ray 2-10 keV luminosity and the obscuring col- 
umn Nh to estimate the nuclear contri bution to the total light 
in each observed band according to the ISilva et all (l2004h tem- 
plates, we subtracted it from the observed photometry, and fitted 
the residual fluxes with the pure stellar library. Stellar masses 
inferred with these approaches agree satisfactorily with those 
obtained through the decomposition method when only Type 
2 AGNs are taken into account. The median/mean value of 
{Mdec ~ M2 fl ,„)/M dec equals to -0.004/-0.046, with a dispersion 
(i.e. semi-interquartile range) of 0.045, while the median/mean 
of (M dec - M agn subtraction) I M dec is equal to 0.040/0.023, with a 
dispersion of 0.190 (these numbers refer to the z = 0.5 - 2.5 
redshift range in GOODS-S). However, both these approaches 
(Moftm and Magn subtraction) fail in reproducing the Type 1 SEDs. 
Indeed, nuclear emission in Syl galaxies does affect all bands, 
not only the near- and mid-infrared ones, and it often outshines 
the stellar light, making the subtraction method unsuccessful. 
The decomposition technique that we adopt is instead the only 
way to properly reproduce Type 1 SEDs (see Fig. [2j». 

Finally, we applied the decomposition technique to non- 
AGN galaxies to check that this procedure is not introducing 
any systematics in the stellar mass estimates. The adoption of 
the stellar+Sy2c77im(c™ci:) library on non-AGNs provides basi- 
cally the same stellar masses, the ratio of those estimated with 
and without decomposition being 0.99 + 0.10 (1.00 + 0.08). A 
larger scatter is obtained when using the stellar+Syl library, the 
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(increasing from bottom to top). Red circles and blue triangles show Type 2 and Type 1 AGN hosts, respectively, and the y-axis 
represents host-only colours. Black small dots identify the inactive population. Vertical dashed lines indicate the mass bins used in 
this analysis. 



previous ratio being equal to 0.98 + 0.32, with only a very mild 
tendency for stellar masses to be smaller when the decomposi- 
tion is applied. However, a 2% error in the stellar mass estimate 
is perfectly consistent with its uncertainty, and represents a fur- 
ther confirmation of the reliability of our mass estimate proce- 
dure. However, in order not to increase the uncertainties, we do 
not include the nuclear component when fitting non-AGN galax- 
ies. 

3.1.1. Colour-mass diagram 

In Fig.|4]we show the colour-stellar mass relation for AGN hosts 
(large coloured symbols) and inactive galaxies (small black dots) 
in the three fields, divided into redshift bins. The figure plots the 
rest-frame (U — V) colour of the best-fit stellar template for each 
object. 

As far as the control sample of non-AGNs is considered, 
thanks to the restriction to good photometric quality, the colour- 
mass diagrams in the three fields are well consistent. 

Figure [4] shows that AGNs are p referentially hosted 



in massive galaxies , as e xpected (e.g. Bundv etall 20081: 
lAlonso-Herrero et all 120081: iBrusa et~aD l2009i: ISilverman et all 



l2009bHXue et al.ll2010l:lMainieri et al.ll201 ll) . although they can 
have a tail towards lower masses. We carefully inspected fits to 
AGNs with suspiciously low stellar mass hosts. We found 8 Type 
1 AGNs in COSMOS with M < 10 10 M o and no 8 fim photom- 
etry. In powerful AGNs, MIR photometry is crucial to constrain 
the nuclear component and recover a reliable stellar mass esti- 
mate. We therefore removed these sources from our sample (they 
are not plotted in Fig. [4] and not accounted when computing the 
statistics given in the previous section). 

The host-only colours of Type 1 AGNs appear very blue 
compared to the host stellar mass in about half of the popula- 
tion. Blue colours with respect to the underlying galaxy popula- 
tion in the host galaxie s of bright AGNs were fi rst reported by 
iKauffmann et alj (120031) and iJahnke et al.l d2004l) . and more re- 
cently bv lCardamone et akl (l2010h . They may hint to some phys- 
ical properties, e.g. the existence of young stellar populations 
with low reddening. However, we also note that large uncertain- 
ties are involved in the SED decomposition of these sources and 
it is difficult to interpret the best-fit stellar colours with any con- 
fidence. The key aspect is that the stellar mass distributions of 
Type 1 and 2 are quite similar and, therefore, the inclusion of 
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Type 1 in our final AGN sample does not bias our subsequent 
results. 



3.2. Far-infrared luminosities and star formation rates 

As a direct tracer of the far-infrared emission we concentrate on 
the mean luminosity vL v (60fim) at a rest-frame wavelength of 
60 fim. This choice is a compromise between a wavelength long 
enough to avoid most of the AGN contamination (see discus- 
sion below) and at the same time short enough to be sampled by 
PACS 160 fim observations even at the highest redshifts consid- 
ered in this work. 

To preclude any assumptions about the SED shape, we com- 
puted vL v (60/im) though a log-linear interpolation of PACS 
fluxes, after converti ng them to luminosities following the ap- 
proach of IShao et all (I2010t) . For the GOODS-S dataset, which 
includes 70 fim observations as well, we interpolated between 
the two PACS bands bracketing rest-frame 60 fim. Although the 
70 fim band is only used in low-z objects (z < 0.67), we checked 
that the addition of these observations in GOODS-S does not in- 
troduce any bias in the analysis: when interpolating rest-frame 
60 fim luminosity only from 100 and 160 fim PACS photometry, 
our results remain unchanged, and vL v (60^m) in each mass and 
redshift bin do not differ by more than 50%. 

We also exp lored fitting PACS flu xe s with a typical IR te m- 
plate (e.g. from lCharv & Elbadl200ll or lDale & Heloull2002l li- 
braries) to derive vL v (60^/m). We found differences by a few tens 
per cent (only in 15% of the bins is the difference larger than 
50%), but the global picture is unaffected. Therefore, we adopt 
the simple interpolated vL v (60fim) estimate throughout, making 
no assumptions on the detailed SED shape. 

For sources fully detected (i.e. detected in both PACS bands 
used to interpolate 60 fim luminosity) we used the correspond- 
ing fluxes and computed individual vL' v (60fim). For PACS un- 
detected sources we computed average fluxes by stacking in a 
given mass and redshift interval. We st acked at the X-ray posi- 
tions on PACS residual maps using the iBethermin et al.l d2010t) 
libraries^. The use of residual maps, from which all detected ob- 
jects were removed, avoids contamination by nearby brighter 
sources. PSF photometry was performed on the final stacked 
images. Uncertainties in the stacked fluxes were computed by 
means of a bootstrap procedure. For each PACS band j, we 
then averaged, by weighting with the number of sources, stacked 
fluxes with individual fluxes of partially detected objects (i.e. de- 
tected only in band j but not in the other one). These stacked and 
averaged fluxes in each band are used to get vL s v TACK {6§fim), in 
the same fashion as for the fully detected sources. 

The final 60 fim luminosity in each mass and redshift interval 
is computed by averaging over the linear luminosities of detec- 
tions and non-detections, weighted by the number of sources. 
Only bins with more than 3 sources are used in our analysis. 

Errors in the infrared luminosity were obtained by bootstrap- 
ping. A set of N sources, where N is equal to the number of 
sources per mass and redshift bin, is randomly chosen 100 times 
among detections and non-detections (allowing thus repetitions), 
and a vL v (60|/m) is computed per each iteration. The standard 
deviation of the obtained vL v (60fim) values gives the error in 
the average 60 fim luminosity in each bin. The error bars thus 
account for both measurement errors (uncertainty in the fluxes) 
and the error in the population distribution (abundances of var- 
ious galaxy types in each bin). They do not, however, account 
for cosmic variance. The stellar mass uncertainties are not ac- 



http://www.ias.u-psud.fr/irgalaxies/downloads.php 



counted for either, since they are small compared to the mass 
bin widths. 

This work is based on the assumption that the rest-frame far- 
infrared emission is dominated by the host galaxy and is there- 
fore a proxy of the SF activity. A number of previous studies 
support o ur assumption based on different grounds . Many au- 
thors (e.g. ISchweitzer et al.ll2006t iNetzer et all2007t iLutz et al.l 
l2008h showed a strong correlation between FIR luminosity 
and SFR tracers, such as PAH emission features, both in lo- 
cal and high redshift bright (Lacn/Lfir up to ~ 10) QSOs. 
lHatziminaoglou et al. (1201 Ol) u s ed a S ED decomposition tech- 
nique based on the Fritz et ail ((2006) dusty torus model, and 
showed that a starburst component is always needed to repro- 
duce both Type 1 and Type 2 AGN FIR e mission. Fo r a de- 
taile d discussion see also the introductio ns of ILutz et al.l (1201 Ol) 
and ISerieant & Hatziminaogloul d2009l) . As a further valida- 
tion, we compared PACS colours (F160/F100) of AGNs with 
those of non-AGNs, and found no significant difference be- 
tween the two samples. Sim ilar results were also obtained by 
lHatziminaoglou et al.l d2010h using SPIRE (250, 350 and 500 
fim) colours. These checks, as well as the results of the stud- 
ies cited above, make us confident that the AGN contamina- 
tion in FIR emission is not dominant over the emission of 
dust heat ed by young sta r s, at l east in the majority of sources. 
However, iMullanev et al.l (1201 ll) suggest that the total emission 
at 60 fim may be dominated by nuclear light in at least three of 
their 1 1 AGN dominated galaxies. In Rosario et al. (in prep.), 
we dem onstrate that, if adopting the templates of lMullanev et al.l 
(l20TTt) . the tail of the torus emission can extend to rest-frame 60 
fim in a fraction of the most luminous (L x > I0 44 erg/s) AGNs. 
We will discuss in Sect. 14. II and 15. II if and how this contamina- 
tion can affect our conclusions. 



4. Results 

In this section we address the main topic of this work, i.e. the 
comparison of average star formation properties of AGN host 
galaxies with inactive galaxies of similar mass. In Fig.|5]we plot 
average vL v (60fim) versus stellar mass for AGN hosts (coloured 
solid lines and circles) and mass-matched non-AGNs (black dot- 
ted lines and open boxes) in the different redshift bins for the 
GOODS sample (South+North, upper panels) and the COSMOS 
one (lower panels). Masses represent the median values in each 
mass and redshift interval. Error bars on the stellar masses show 
the 80% confidence intervals. 

Given the different depth and sky areas, the GOODS and 
COSMOS fields can be regarded as somewhat complementary 
samples (see Fig. [1). However, given the large difference in 
the IR, X-ray and optical depth between COSMOS and the 
GOODS fields, combining the catalogues would introduce bi- 
ases which can not be easily taken into account. For these rea- 
sons we decided to keep the COSMOS analysis separate and to 
combine GOODS-S and GOODS-N, which have similar prop- 
erties and show similar results if considered separately, to in- 
crease the statistics in the deep fields. In order not to stack 
on different PACS residual images with different noise levels, 
vL v (60fim) was independently computed in each redshift and 
mass bin for each of the GOODS fields. We then computed the 
mean GOODS vL v (60fim) for each bin by averaging the val- 
ues obtained separately in the two fields, weighting them by 
the corresponding errors. Similarly, the typical mass associated 
to a given vL v (60fim) was computed by averaging the individ- 
ual median masses weighted with the errors in the individual 
vL v (6Qfim). The inter-percentile range of the stellar mass dis- 
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Fig. 5. Average vL v (60/im), proportional to the SFR, versus stellar mass in the different redshift bins, increasing from left to right, 
for the GOODS fields (upper panels) and COSMOS (lower panels). Coloured solid lines/solid circles and black dotted lines/open 
boxes show AGN hosts and inactive galaxies, respectively. Stellar mass values reflect the median in each mass and redshift interval, 
and their error bars show the 80% confidence intervals. Error bars on vL v (6Q ^m) were computed by bootstrapping (see text). Grey 
dashed lines represent the main sequenc es as inferred bv ISantini et all d2009h in r edsh ift bins simi l ar to the ones used in this work, 
while dot-dashed ones are the results of iNoeske et al.l d2007l) . lElbaz et all d2007h and lDaddi et al.l (l2007bi) respectively at z ~ 0.5, 
Z ~ 1 and z ~ 2. At all masses and redshifts, with the exception of the highest redshift interval in the GOODS fields, AGN hosts 
show an enhanced FIR emission with respect to inactive galaxies of similar stellar mass. In contrast, the locus of AGN hosts is 
broadly consistent with the main sequence of (only) star forming galaxies. 



tribution at 80% was computed on the full combined GOODS 
sample in each bin as an estimate of mass uncertainty. 

Each field is affected by incompleteness at low masses in 
a different way. In the mass range studied (10 9 - 10 12 M o ), our 
control sample lacks sources in the GOODS-S field at the high- 
est redshift bin, which is only complete above ~ 2.5 ■ 10 9 M Q , 
and in COSMOS, where the limit is set around ~ 4 • 1O 9 M at 
0.8 < z < 1.5 and ~ 5.5 ■ 1O 9 M at 1.5 < z < 2.5. However, mass 
incompleteness is not a serious issue for this kind of work for 
two reasons. Firstly, we only deal with mass-matched samples, 
which, given the AGN mass distribution, only include extremely 
few sources (only 6 AGNs in COSMOS) below the mass incom- 
pleteness limits. Secondly, the aim of this study is not measuring 
the absolute SF activity, but rather comparing SF properties be- 
tween X-ray selected AGN hosts and inactive galaxies of similar 
mass, and mass incompleteness can be ignored as long as it af- 
fects the two samples in a similar way: since in AGN hosts the 
rest-frame /if -band, which is broadly proportional to the stellar 



mass, is generally dominated by stellar emission, there is no rea- 
son to believe that this is not the case. 

Using the 60 fim rest-frame luminosity as a proxy of the star 
formation rate, a vL v (60jum)-stellar mass diagram is analogous 
to one plotting SFR-stellar mass. From Fig. [5] it is possible to 
observe the global increase of the far-infrared luminosity, and 
hence SFR, with redshift as we ll as with stellar mass, in both ac- 
tive and inactive galaxies (e.g . lSantinT et al. 20091: iDamen et all 
12009b iRodighiero et alJl2O10[ and references therein). We also 
plot the main sequence relations (grey curves) inferred by pre- 
vious stud ies in similar redsh ift intervals as the ones used in 
this work: ISantini et al] (120091) at 7 = 0.6 - 1.0, z = 1.0 - 1.5 
and z = 1 . 5 - 2 .5, |No eskeetal] d2007| ) at z = 0.2 - 0.7 
lElbaz et al] d2007l) at z = 0.8 - 1.2 and iDaddi etal] d2007bl) 
at z = 1.4 - 2.5. For this purpose, we first converted the SFR 
into an estimate of the total infrared luminosity L/fl(8-iooo/im) by 
dividing by 1.8 10" 10 (following ISantini et al] [20091) and then 
converted the latter into vL v (60fim) by linearly fitting the values 
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of £/j» 8-ioooum) and vL v (60yum) predicted by the ICharv & Elbazl 
d200lh library. 

Our control sample of inactive galaxies shows general agree- 
ment with the locus of the main sequence with the exception of 
the highest masses. The slopes for our sample indeed appear flat- 
ter than the main sequence. This is explained by the fact that the 
main sequence relations are computed for purely star forming 
galaxies, while our mass-matched reference sample is made by 
a combination of star forming and quiescent galaxies. The latter, 
more abundant at high masses and low redshifts, are responsible 
for the bending of the relation. 

One may also notice that the locus of inactive galaxies (black 
curves) in both the average GOODS field and the COSMOS field 
are quite similar, especially at intermediate masses (between 
10 10 - 10 n M Q ) where we have the best statistics. At higher 
masses, cosmic variance can lead to the small differences we ob- 
serve. At low masses, incompleteness in our samples can play a 
role. However, despite these effects, the trends for inactive galax- 
ies between fields are generally consistent, despite their very dif- 
ferent PACS depths. This may be taken as evidence for the va- 
lidity of the procedure we employ to derive vL,,(60/mi). 

The locus occupied by AGN hosts is also roughly consis- 
tent with the main sequence, with the exception of the lower 
mass bins which show a larger SFR. They could be pushed up 
by a higher fraction of strongly star forming galaxies at low 
masses. Broad consistency with t he main sequence wa s also 
reported by the previous studies of iMainieri et al. l d20Tlh . who 
measured comparable specific SFRs in obscure d quasars and in 
an IRA C 3.6 yum selected control sample, and iMullanev et all 
d2012l) . who found that the specific SFRs in AGN hosts with 
Lx = 10 42 - I0 44 erg/s are only marginally lower (by 20%) than 
those of main sequence galaxies at 0.5 < z < 3.0. 

The main target of this study is however the difference be- 
tween AGN hosts and a mass-matched control sample of inac- 
tive galaxies, which, being optically-selected, includes both the 
star forming and the quiescent populations. A comparison sam- 
ple including both star forming and quiescent galaxies is mo- 
tivated by, e.g., the wide spread of AGN loci in optical color- 
magnitude diagrams. The measured vL v (60/mi) in each redshift 
and mass bin and in each field are reported in Tab. [2] Both from 
the latter table and from Fig. [5] it is clear that at all redshifts and 
masses AGN hosts sit above non-AGNs, with the exception of 
the GOODS high redshift interval, where the two samples show 
consistent FIR emission. The difference in the FIR emission be- 
tween the two populations is larger in COSMOS, and more mod- 
est in GOODS. We also note that the vL v (60/mi) enhancement 
at low redshift seems to be stronger than at high redshift. This 
can be partly explained by the slightly larger (0.1 - 0.2) aver- 
age redshift of the inactive control sample compared to AGNs at 
z = 1.5-2.5 (especially in GOODS-N, see Sect.l2~4ii, and makes 
the global offset that we observe between the two populations 
even m ore robust. I t is al so in qualitative agreement with the re- 
sults of IXue et al.l (1201 Oh . who report a larger difference in the 
SFR between AGN hosts and mass-matched non-AGN galaxies 
in GOODS at z - - 1 (~ 0.3 - 0.4 dex) than at z - 1-3 
(~ 0.1 -0.2 dex). We report in Tab.[3](second and fifth columns) 
the average vL v (60fim) offsets^ in each redshift bin and in the 
total redshift range. The FIR emission in AGN hosts is modestly 
enhanced in the GOODS fields. The average offset is equal to 



0.32 + 0.03 dex (the average offsets in the two GOODS fields 
considered separately is equal to 0.30 + 0.03 dex in GOODS-S 
and 0.17 + 0.03 dex in GOODS-t-fl)- In COSMOS the level of 
enhancement is larger, with an average offset of 0.42 + 0.04 dex. 
We will discuss in Sect. 15.11 how these different results found 
for the GOODS vs COSMOS fields are driven by the different 
luminosities of the AGN samples (Fig.Q]). 

Enhanced star formation in high z AGN hosts may at first 
glance be surprising given the tendency of these hosts to pref- 
erentially reside in the so- called green valley i n optical colour- 
magnitude diagram s (e.g. iNandra et alj|2007t : ICoil et all 120091 : 
iHickox et al.l 120091) . It has been shown, however, that this ten- 
dency is strongly influenced by varying mass-to-light ratios 
and disappears when study i ng colours of m ass-selected sam- 
ples dSilverman et al.|[2009ab IXue et ai1l2010l) or even reverts t o 
a preference for blue cloud/green valley hosts dAird et alj|2012l) . 

4.1. A closer look at SFR distributions 

So far we have discussed average SF properties of the AGN 
host population compared to a control sample of inactive galax- 
ies of similar mass. However, averages are insensitive to the 
detailed distribution of SFRs of the two samples, especially 
if they are not normally dist ributed, as is known to be the 
case for inactive ga l axies ( e.g. Bell et al] 2004t ISalimbeni et al.l 
l2008b ISantini et all I20091: lllbert et al.l 120101) and AGN hosts 



(e.g.lNandra et alj2007M Silverman et al.l2008t[Cardamone et al.l 
I2010ir To gain a handle on the underlying SFR distributions 
of the two samples, we plot in Fig. [6] vL v (60/jm) versus stellar 
mass for individual PACS detections (small symbols) as well as 
the stacked values for non detections (large symbols), both for 
AGNs (coloured solid circles) and for the control sample (black 
open boxes). 

We distinguish between the low-Lx subsample (on the top), 
where the FIR luminosity can be safely regarded as a proxy of 
the SFR (see discussion in Sect. 13. 2t . and the total sample (on the 
bottom), where nuclear light may affect the vL v (6Q/-im) estimate 
in some of the most luminous AGNs. In order to split the sample 
by X-ray luminosity, we chose redshift-dependent Lx thresholds, 
equal to log Lx[erg s~'] =43.3,43.8 and 44.0 respectively at z ~ 
0.65, 1.15 and 2, shown as red lines in Fig.[T](for a justification, 
please see Sect. 15. U . As far as the low-Lx subsample is consid- 
ered, a new set of control galaxies was extracted from the parent 
inactive sample to match the mass distribution of low luminosity 
AGN hosts. 

First of all, it must be noted that, in the deeper GOODS 
fields, individual PACS detections are able to probe over a good 
fraction of the main sequence of star forming galaxies, at least 
at higher masses. In shallower fields like COSMOS only the up- 
per envelope is traced by individual detections, and we have to 
rely on stacks to investigate the lower envelope (as well as the 
quiescent population). 

Among PACS detected sources, no significant offset in FIR 
luminosity is observed between AGN hosts and inactive galax- 
ies. This can be seen from the inset histograms in Fig. [6] al- 
though the low number of detected galaxies in some bins does 
not allow a proper comparison. According to a Kolmogorov 
Smirnov test performed on both the total sample and the low- 
Lx subsample, the log(vL v (60yum)) distributions for PACS de- 



3 We computed linear averages of the differences between the loga- 
rithmic FIR luminosities in the individual bins; uncertainties in the av- 
erage were computed by summing in quadrature the logarithmic error 
bars; upper limits were not considered in the average computation. 



4 It is unsurprising that the two offsets considered separately are 
lower than the average offset. Indeed upper limits, more easily in place 
when dealing with only one field and the number of sources is lower, 
are not considered in the average computation. 
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Table 2. Average log(vL v (60yum)[erg s -1 ]) values in each redshift (rows) and mass (columns) bin in the different fields. 
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Table 3. Average FIR emission offsets between AGNs and non-AGNs over all mass bins. 



Average offsets 



GOODS COSMOS 





z 


Total sample 


Low-Lx High-Lx 


Total sample 


how-Lx 


High-L x 


0.5 


-0.8 


0.41 ± 0.04 


0.29 ±0.06 


0.44 ± 0.06 


0.21 ± 0.05 


0.81 ±0.07 


0.8 


- 1.5 


0.60 ± 0.06 


0.56 ± 0.06 


0.44 ±0.10 


0.50 ± 0.09 


0.45 ± 0.04 


1.5 


-2.5 


-0.06 ± 0.05 


-0.01 ±0.06 


0.37 ± 0.04 


0.02 ± 0.07 


0.42 ± 0.04 


0.5 


-2.5 


0.32 ± 0.03 


0.28 ±0.10 


0.42 ± 0.04 


0.24 ± 0.04 


0.56 ± 0.03 



Notes. The different rows correspond to different redshift bins, while the columns correspond to the entire sample, the low- and the high-L^ 
luminosity ones. Average logarithmic offsets were computed as linear averages of the differences between the logarithmic FIR luminosities. 
Uncertainties were computed by summing in quadrature the logarithmic error bars on vL v (60jim). Upper limits were not considered in the average 
computation. Bins with less than three sources were discarded. 

Table 4. Detailed comparison between AGNs and non-AGNs by separately considering PACS detected and undetected sources. 









GOODS 






COSMOS 






z 


(1) 


(2) 


(3) 


(1) 


(2) 


(3) 










Total samp] 


c 






0.5 
0.8 
1.5 


-0.8 
- 1.5 
-2.5 


77% 
1% 
82% 


0.28 ± 0.04 
0.38 ± 0.06 
-0.13 ±0.04 


3.87 ± 1.00 
2.02 ± 0.99 
1.76 ± 1.00 


64% 
18% 
36% 


0.43 ± 0.06 
0.22 ± 0.06 
0.28 ± 0.04 


2.27 ± 1.00 

3.28 ± 1.02 
10.17 ± 1.30 


Low-Lx subsample 


0.5 
0.8 
1.5 


-0.8 
- 1.5 
-2.5 


76% 
1% 
59% 


0.22 ± 0.04 
0.39 ± 0.05 
-0.14 ±0.05 


7.26 ± 1.29 
1.71 ±0.98 
1.44 ± 1.39 


99% 
83% 
31% 


0.23 ± 0.05 
0.17 ±0.08 
-0.08 ± 0.08 


1.90 ± 1.03 
3.78 ± 1.06 
4.38 ± 1.41 



Notes. Columns: (1) probability, according to a Kolmogorov Smirnov test, that the log(vL, / (60/jm)) distributions of PACS detected AGNs and 
non-AGNs are consistent with describing the same population; (2) average offset between the stacked log(vL v (60/jm)) for PACS undetected AGNs 
and non-AGN; (3) ratio of PACS detection rates of AGNs and non-AGNs. The upper table refers to the low luminosity subsample (upper panel of 
Fig.|6j, whereas the lower table represents the total sample (lower panel of Fig.|6). 



tected AGNs and non-AGNs are consistent with describing the 
same population of galaxies in most of the bins (the probabil- 
ity is reported in column (1) of Tab. |4}. In the GOODS inter- 
mediate redshift interval, these probabilities are biased by three 
or four outliers: when removing such sources, the probabilities 



that the two distributions describe the same populations turn out 
to be 10% and 3%, respectively for the low-L^ subsample and 
the total sample. In the low-L^ subsample, unbiased by possi- 
ble AGN contamination to the FIR, the distribution of PACS de- 
tected AGN hosts and non-AGNs are therefore consistent at all 
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redshifts and in both fields. The slightly lower probabilities mea- 
sured in the COSMOS total sample can be attributed to either an 
intrinsically larger SFR in high-Lx PACS detected AGN hosts or 
nuclear contamination to vL,,(60 / um) in the brightest AGNs. 

An offset between AGN hosts and non-AGNs is instead ob- 
served among the stacked non detections, including weak star 
forming as well as quiescent galaxies (see column (2) of Tab. 
0). The averag^l offsets at all masses and redshifts between the 
stacked vL v (60^m) of AGNs and the control sample for low/all 
X-ray luminosity is 0.16 + 0.03/0.18 + 0.03 dex in GOODS and 
0.11 +0.04/0.31 +0.03 dex in COSMOS. However, these offsets 
are not large enough to completely explain the SF enhancement 
observed in Fig. This seems to exclude that a significant role 
may be played by a larger fraction of quiescent galaxies in the 
control sample, which would imply larger offsets among unde- 
tected sources. 

Since averages are computed by weighting with the num- 
ber of sources, it immediately follows that PACS detection rates 
must favour AGN hosts. To verify this statement, we plotted in 
Fig.|7]the fraction of PACS detected sources among AGN hosts 
(coloured solid lines) and inactive galaxies (black dotted lines) in 
each of our redshift and mass bins. As before, we required PACS 
detections in both bands used for vL v (60/im) interpolation in or- 
der for an object to be 'detected'. AGN hosts of low/all X-ray 
luminosities are on averag eS 3.5 +0.7/ 2.6 +0.6 times more likely 
in GOODS and 3.4 + 0.7/5.2 + 0.6 times more likely in COSMOS 
to be detected by PACS than the control galaxies (the ratios in 
each redshift bin are reported in column (3) of Tab.|4|i. The PACS 
detection rates in each redshift and mass bin for the total sam- 
ple (shown in the right panel of Fig. 13 can be found in Tab. 
Q] We also checked this results against possible AGN contami- 
nation at the shortest PACS bands. The excess in the detection 
rates of AGN hosts compared to non-AGNs is the same within 
the errors if detection is restricted only to the PACS 160^111 or 
PACSlOOpm bands, in agreement with the fact that AGN hosts 
and inactive galaxies have consistent F160/F100 colours (see 
Sect. 13.2b . This further strengthens the assumption that AGN 
contamination does not strongly influence the vL v (60pm) deriva- 
tion, even for the total sample, and validates the interpretation of 
vL v (60^m) as a SFR diagnostic. 

We can conclude that the SF enhancement that we observe in 
Fig. [5] is on average not attributable to larger SFRs in rare very 
extreme starburst events in AGN hosts than in inactive galax- 
ies, which would shift the SFR distribution to higher values. It 
is instead due to a combination of modestly brighter stacked 
FIR emission in faint FIR sources (i.e. undetected by PACS) 
and higher PACS detection rate (i.e. larger fraction of highly star 
forming galaxies) among AGN hosts. By means of a simple cal- 
culation described in Appendix A, the higher PACS detection 
rate turned out to be the dominant effect in the deep GOODS 
fields. Indeed, for a given detection rate, the linear average is 
highly dominated by the brightest sources, with PACS unde- 
tected sources only contributing in a minor fashion to the aver- 
age offset. In other words, the distribution of SFRs among AGN 
hosts is not grossly different from norm al star forming galaxies , 
consistent with the principal result from lMullanev et all (|2012|) . 



5 Average offsets were computed in the same way as for the global 
(detections+non-detections) averages, and upper limits were not con- 
sidered. 

6 We computed linear averages from the individual bins of the ratios 
between the detection rates of the two populations; uncertainties in the 
average were computed by summing in quadrature the binomial uncer- 
tainties in the individual detection rates; upper limits were not consid- 
ered. 



Instead, on average, star forming galaxies are more likely to 
host AGNs at a given stellar mass. However, in the shallower 
COSMOS field, given the larger fraction of PACS undetected 
sources, the brighter FIR stacked emission and the higher PACS 
detection rate among AGNs can give comparable contributions 
to the average vL v (60/im) offset between the two populations. 

5. Discussion 

From the results presented in the previous section, it is clear that 
AGN hosts have enhanced FIR emission with respect to a mass- 
matched sample of inactive galaxies including both star forming 
and quiescent populations. This is observed in all fields and at all 
masses and redshifts except the highest redshift bin in GOODS. 
The level of SF enhancement between the two populations de- 
pends on the fields, with COSMOS showing a larger offset than 
GOODS. As shown in Fig. [Q COSMOS probes to higher X- 
ray luminosities compared to the GOODS fields, while missing 
most low-L x AGNs. In some ways, COSMOS and GOODS are 
complementary AGN samples. This hints to the possibility that 
the level of SF enhancement is related to the X-ray emission 
from the AGN. One interpretation, following the scenario sug- 
gested by our pre vious resul t s infer red with SDP Herschel data 
and presented in IShao et al ] (I20101) . is that the different levels 
of SF enhancement observed in the two fields reflect different 
modes of AGN growth. However, before exploring this possibil- 
ity, we mention that another viable interpretation of the larger 
SF activity observed in AGN hosts is in terms of positive feed- 
back; AGN outflows, such as winds and radio jets, compress 
and collapse interstellar clouds, thus inducing SF, as suggested 
by a numb er of previo us wo rks (e.g. iBegelman & Cioffill l989; 
[Sil5l2005l iFeain et al.ll2007l ISilk & Normanll2009t lElbaz et alj 
120091) . Radio-based studies of this population may help to un- 
cover the role of positive feedback. 



5.1. Two modes of AGN evolution? 

To verify whether the different vL v (60/iin) offsets observed in 
GOODS and COSMOS are a consequence of the different AGN 
luminosity regimes mostly probed by each field, and hence pos- 
sibly a consequence of different AGN evolutionary modes, we 
repeated the analysis described in Sect. [4] by dividing the sam- 
ples into high- and low-Lx subsamples by using the redshift- 
dependent Lx thresholds reported in Sect. I4.ll As done for the 
complete sample and as mentioned above, we extracted mass- 
matched subsamples from the total control dataset to avoid any 
bias due to the different stellar mass distribution of AGN hosts 
and inactive galaxies. 

Another motivation for splitting the dataset based on the X- 
ray luminosity is to arrive at samples free of possible AGN con- 
tamination to the FIR emission: in the low-Lx host galaxies, the 
vL v (60yum) enhancement can be safely regarded as a SFR en- 
hancement (see discussion in Sect. 13. 2\ . 

The results are presented in Fig. [8] and the average 
vL v (60yum) offsets are summarized in Tab. [3] Once the same 
AGN luminosities are sampled, GOODS and COSMOS fields 
provide consistent results. The GOODS fields are almost entirely 
(~ 88%) made of galaxies fainter than the luminosity thresh- 
olds adopted to split the sample. Consequently, the average 
vL v (60yum) offset observed in the low-Lx subsample (0.28+0.03 
dex, top left panels) does not differ significantly from that ob- 
served in the total sample. The high-Lx subsample does not of- 
fer enough statistics to extract meaningful numbers, owing to 
the paucity of luminous sources in the small GOODS fields. In 
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COSMOS, the average FIR offset measured in the high-Lx sub- 
sample (bottom right panels) increases to 0.56 + 0.03 dex, while 
at lower X-ray luminosities (bottom left panels) it is reduced to 
0.24 + 0.04 dex, consistent with the FIR enhancement observed 
in GOODS at similar AGN luminosities. 

The consistency between GOODS and COSMOS results, 
when the same AGN population is considered, is a further con- 
firmation that uncertainties involving the host stellar mass esti- 
mate in Type 1 AGNs are not introducing serious biases in the 
analysis. Indeed, while the GOODS low-L v subsample is almost 
completely (~ 90%) made of Type 2 AGNs, the fraction of Type 
1 AGNs in the COSMOS low-Lx one is ~ 30%. It also ensures 
that the analysis in the COSMOS field is insignificantly affected 
by the /-band selection applied to the sample (Sect. l2~4l i. 

ISilverman etaD (l2009bl) used [O II] (A = 3727 A) as a star 
formation indicator to look for enhanced star formation in z < 1 
AGN hosts. They found a weak but not significant enhancement 
for hosts of moderate luminosity (10 42 < Lx[erg s _1 ] < 10 43 7 ) 
AGNs, while indication for enhanced SF if considering hosts of 
more luminous AG Ns. Intricacies of [O II] as an AGN host star 
formation indicator dKim et al.ll2006h likely reconcile this with 
our findings. 

In order to understand whether the different vL v (60/im) en- 
hancement that we observe in the two subsamples is indeed 
due to the different accretion regimes, a few effects must be 
taken into account. Although a morphological analysis is beyond 
the scope of the present study, we note that the morphological 
mix in the low-Lx a nd high-Lx AGNs subsamples may in gen- 
eral be different (e.g.lHutchings et al.ll2002Ll2009t IVeilleux et all 
l2009atlMainieri et al.l201 UlKocevski et al.l2012l) . The same can 
in principle be true also for the mass-matched control galaxies, 
even though the mass distribution of low-Lx and high-Lx AGN 
host galaxies are very similar, with the latter having an average 
stellar mass ~ 0.14 + 0.03 dex higher than their lower luminos- 
ity counterparts (the comparison has been done in the COSMOS 
field, where we have the better statistics). Since the combina- 
tion of these effects is hard to account for, we compared the FIR 
emission among low-Lx vs high-Lx AGN host galaxies, and, 
similarly, among the inactive galaxies mass-matched to low-Lx 
AGN hosts vs those mass-matched to high-Lx AGN hosts. As 
can also be seen from Fig.[8j the two control samples are totally 
consistent, while an offset is observed among low-Lx and high- 
Lx AGN hosts, the latter showing the larger FIR emission. The 
different behaviour of the two subsamples can therefore be asso- 
ciated with differences in the SF activity among the host galaxies 
of faint and bright AGNs. 

The Lx-dependent FIR enhancement can be interpreted in 
terms of the AGN evoluti onary picture out lined by the previous 
step of the present work dShao et al.ll2010h and described in the 
introduction. There, we used Herschel SDP data from GOODS- 
N to study the relation between host star formation and bolo- 
metric AGN luminosity at different redshifts and we showed that 
AGN growth can be explained by a combination of two paths of 
evolution. Very luminous AGNs are closely coupled to their host 
galaxy's growth, through a process that leads to a correlation be- 
tween AGN luminosity and SF activity in the host determining a 
tight co-evolutionary phase. A good possibility are gas-rich ma- 
jor mergers: these transport to the centre large amounts of gas 
which trigger both AGN fuelling and intense SF bursts. On the 
other side, at lower luminosities AGN hosts experience secu- 
lar evolution, in the sense that instantaneous nuclear accretion 
and total secular, non-merger driven SF are not strongly corre- 
lated and expected to evolve independently. In Rosario et al. (in 
prep.), which represents our third step in the study of the inter- 



play between AGN accretion and host's star formation, we ex- 
tend this analysis by including more Herschel fields and make 
use of a larger dataset which allows sampling of both low and 
high luminosity AGNs. We find similar results wit h respect to 
our pr evious study, although, as already suggested bv lShao et al.l 
(I2010I) . the transition between the two modes of AGN growth 
seems to be shallower at high redshift. This transition is mea- 
sured around Lagn ~ 10 44 5 ~ 45 5 erg s _1 which correspond^ 
to Lx ~ 10 43 - 3-44 erg s _1 . These redshift-dependent luminosity 
thresholds were used to split the sample and are shown as red 
lines in Fig.Q] 

The two-mode evolutionary scenario above explains the 
larger vL v (60/im) enhancement (with respect to mass-matched 
inactive galaxies) observed in high-Lx AGN hosts. However, we 
remind that contamination from the AGN could affect the FIR 
luminosity estimate in a fraction of the brightest (Lx > 10 44 erg 
s~') hosts (see Rosario et al. in prep, for further details). The 
aforementioned shallower transition between the two growth 
modes observed at high redshift also explains the larger SF en- 
hancement measured at low z than at high z. This behaviour is 
expected for two reasons: first, given the lower average SFR at a 
given mass, we are more sensitive to merger-induced SF bursts; 
second, larger errors at high z and possible gas feeding through 
secular processes (besides major mergers) can wash out the ob- 
served SF offset. 

On the other hand, the IShao et all d2010l) picture does not 
predict any vL v (60^/m) enhancement, i.e. SF enhancement, in 
low-Lx AGN hosts. In this view, the offset in low-Lx is there- 
fore rather surprising. We discuss possible interpretations in the 
next section. 

There is another possibility to explain the difference in the 
FIR excess observed with respect to inactive galaxies between 
hosts of low-Lx and high-Lx AGNs, without assuming the ex- 
istence of separate modes of evolution (i.e. merger co-evolution 
vs secular). AGNs may be more likely to exist in galaxies with 
higher gas fractions at a given stellar mass, only because cold 
gas in these galaxies has a higher probability of reaching the nu- 
cleus. Indeed, in order to feed the AGN, a substantial amount 
of gas is required to inflow to the centre, plausibly from large 
scales. Low-Lx would require smaller gas reservoir than brighter 
AGNs, and hence would be characterized by less star formation 
activity on large scales. Under this hypothesis, we would ex- 
pect to see a correlation between AGN and FIR luminosity at all 
scales, while there is no evide nce for such trend at log(Lx[erg 
s -1 ]) < 43.5 dShao et al.ll2010l Rosario et al. in prep.). For this 
reason, we d ecided to in t erpret our data in the view of the picture 
presented bv lShao et al.l d2010l) . However, deeper data and larger 
statistics is needed to definitely accept or exclude this hypothe- 
sis. 



5.2. SF enhancement in low-L x AGN hosts: 

non-synchronous accretion/SF histories or break in the 
dichotomy? 

While the FIR enhancement in bright AGN hosts can be ex- 
plained by merger-induced SF bu rsts, the twofold A GN evolu- 
tionary path that we suggested in IShao et al l d20Tol) cannot ex- 
plain any SF enhancement in low-Lx AGN hosts. However, a 
modest but significant offset of ~ 0.26 dex at ~ 3cr confidence 
level is measured with respect to the mass-matched inactive con- 



7 Bolometric AGN luminositi es are converted into in trinsic 2-10 keV 
luminosities by using Eq. (5) of iMaiolino et al.l d2007l) and adopting a 
ratio of 7 between bolometric and 5 100 A luminosity. 
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trol sample. We suggest here some possible reasons for this en- 
hancement. They may be broadly classified as enhancements in 
a subset of AGN hosts, due to gas-rich major galaxy mergers, or 
enhancements in most AGN hosts, due to secularly driven rela- 
tionships between accretion and star formation. 

The first group of possibilities implies that in objects fol- 
lowing a 'merger' co-evolution, the instantaneous AGN accre- 
tion and star formation rate are not well synchronized, even if 
causally linked. This will lead to comp lex evolutionary t racks in 
a Lpir/Lacn diagram (such as Fig 6 of Sh aoet al.ll2.0ldh . and to 
merger-enhanced SFRs being present in objects with low AGN 
luminosity. Specifically, it is plausible to discuss the case of de- 
layed AGN feeding with respect to the onset of the SF episodes 
induced by the gas inflow during major intera ctions. This cir- 
cums t ance is suggested by s everal models (e.g. iDi Matteo et ail 
120051: iHopkins et all 120061: iNetzed l2009h . which assume that 
AGNs are fuelled only when substantial amounts of gas can 
reach the central black holes. It is in line with the classical local 
ULIRG-to-QSO transitio n scenario (ISanders et al.lll988bl fl989t 
ISanders&Mirabelll996l) . and would concern a small number of 
strongly SF galaxie s. It is also supported b y simulations such as 
those presented by lHavward" et al.) d201lh. which are based on 
the observational evidences provided by Alexander et al.l (l2008h 
and show that accretion/SF histories in SMGs are not strictly 
synchronized. iNetzerl (120091) supports this scenario by arguing 
the existence of sources, such as ULIRGs and SMGs, character- 
ized bj^JrugeS>FR^_and modest AGN contribution (see also Fig. 
6 of iLutz et alJl20l6T) . According to the evoluti onary sequenc e 
for the starburst-AGN co-evolution suggested bv lNetzerl (120091) . 
the SF episode starts first, and only after a certain time lag (At$ f) 
does enough cold gas reach the centre and the black hole accre- 
tion begin. Simulations prefer a long Ats f compared to the rise 
time of the AGN phase (A^caO- Such a model is qualitatively 
compatible with the SF offset we measure only if AGN duty cy- 
cles are long with respect to the SF episode (i.e., AtAGN 5= Af^f ). 
This will lead to a number of merger-induced starbursts with low 
levels of accretion, which would raise the mean SFR of low lu- 
minosity AGN hosts compared to inactive galaxies. If, on the 
other hand, Af^c/v << Ats f, a number of the strongest star form- 
ing galaxies would be missing from the AGN sample, leading 
to a higher average SFR for non-AGNs, at variance with what 
we find. Alternatively, one may relax the requirement that bright 
AGNs occur only after most of the star formation has progressed 
in mergers, and instead postulate that the general synchroniza- 
tion between the two phases is rather weak. While this can quali- 
tatively explain the observed mean offset, it is not well-motivated 
by most models of AGN fuelling by tidal torques. 

There is a further explanation for the low but significant SF 
enhancement observed at low X-ray luminosities. Alternative 
AGN feeding mechanisms, at work during secular evolution and 
not invoking major mergers, may lead to a weaker physical con- 
nection, but still a certain level of correlation, between instan- 
taneous accretion rate and total star formation. This will make 
the dichotomy between secular evolution and co-evolution less 
sharp than outlined in Sect. 15.11 The nature of the relationship is 
such that the level of accretion is not intimately correlated with 
the level of star formation; it may be causal or may be mediated 
through other processes. For example, one may imagine star for- 
mation in the circum-nuclear environment of SMBHs can dis- 
turb gas by supernovae or stellar winds and lead to inflow to the 
nuclear regions. In this picture, only star formation within the 
inner kpc or so may be related to AGN activity. Alternatively, 
processes that drive star formation in inactive galaxies, such as 
fresh gas infall, bar instabilities, stellar feedback, galaxy harass- 



ment or satellite accretion, may be able to inspire some inflow 
to the SMBH, but to a degree that is not closely related to the 
level of star formation at large scales. An example of these pro- 
cesses, which applies to z > 2 host s, is gravitational insta bility 
in c lumpy disks. The recen t work of iBournaud et al.l d201 lb (see 
also lDi Matteo et al.l20lTl for bright OSOs atz > 5) suggests that 
both AGN accretion and bulge building, often characterized by 
giant star forming clumps, can be associated with violent grav- 
itational instability in high redshift disk galaxies. This can be 
generated by cold strea ms and high gas fra ctions, and it includes 
clumps migration (e.g. lGenze l et al. 2008). A prediction of their 
simulation, which is supported by observational evidences (see 
below), is that high-z moderate luminosity (Lx ~ 10 42 ~ 43 erg Is) 
AGNs are hosted by star forming clumpy disk galaxies with 
no signs of major mergers. Finally, dry mergers, with low gas 
fractions, could also in principle explain the reduced SF en- 
hancement in low-Lx AGN hosts. How ever, most low luminos- 
ity AGNs are in disks at these redshifts dSchawinski et al.ll201 ll 
iKocevski et al.1 120T2I) . so it is unlikely that they underwent a 
(dry) merger recently. 

According to the first class of explanations, a subset of low- 
L x AGNs are associated with strong merger-induced starbursts, 
while the rest are in normal star forming galaxies. Therefore, we 
would expect low-L^ AGNs to be over-represented at the upper 
envelope of the star forming main sequence, where most star- 
bursting systems lie. Our analysis in Sect. 14.11 suggests, how- 
ever, that bright FIR AGN hosts and inactive galaxies have 
rather similar distributions in the SFR-stellar mass plane. Thus, a 
merger-driven scenario for the offset we measure is not strongly 
supported by our observations. The lack of an excess of ma- 
jor mergers among low luminosity AGNs is consistent with 
several observational studies, according to which moderate lu- 
minosity AGNs have similar morphologies to in active galax- 
ies, not consistent with d isturbed morphologie s dGrogin et al.l 
120051: iGabor et al.1 120091: ICisternas et ail 1201 ll). and are fre- 
quently found in disk s at high redshifts dSchawinski et al.ll20lTt 
IKocevski et al.1 120121) . The latter property is consistent with 
AGNs being on av erage hosted by main sequence galaxies (e.g. 
IWuvts et al.l 201 lb . The second class of scenarios applies to the 
majority of normal star forming galaxies and predicts that AGN 
hosts have similar star formation properties as inactive galaxies, 
but that AGNs are a bit more likely to be found in hosts where 
some degree of recent star formation is in effect. This is qualita- 
tively consistent with our study of the SFR distributions, as well 
as the fairly constant SFR that we me asured in low-L^ AGNs 
across a range of Lx dShao et al.ll2010l) . However, further statis- 
tics, deeper FIR and multiwavelength observations and careful 
morphological studies are needed in order to definitely distin- 
guish between the two hypotheses, understand the importance 
of mergers as AGN hosts, as well as explore the processes of 
secular fuelling in AGNs at these redshifts. 

6. Summary 

We used Herschel PACS observations carried out in the two 
GOODS fields and in COSMOS, in combination with ancillary 
optical-to-MIR and X-ray data, to study the star forming proper- 
ties of AGN hosts as a function of stellar mass and redshift with 
respect to inactive galaxies. A fitting procedure to estimate the 
stellar masses of galaxies was developed to disentangle stellar 
and nuclear emission, in order to measure stellar masses of host 
galaxies unbiased by the AGN light. The FIR emission, which 
can be considered as a proxy of the SFR in all galaxies except 
in a fraction of the brightest AGNs, was derived by interpolating 



16 



P. Santini et al.: Enhanced SFR in AGN hosts from PEP-Herschel observations 



PACS fluxes to derive a rest-frame 60 fim luminosity, without as- 
suming a particular FIR SED shape. The FIR emission of AGN 
hosts was compared to that of non-AGN control galaxies with 
similar masses and at similar redshifts including both star form- 
ing and quiescent populations. 

We found that the FIR emission, and hence the SF activ- 
ity, in AGN hosts is enhanced with respect to that of mass- 
matched inactive galaxies. We checked our analysis against pos- 
sible sources of bias (PACS source extraction techniques, dif- 
ferent spectroscopic fractions among the various samples, large 
uncertainties in the host stellar mass estimates in Type 1 AGNs, 
different methods to measure the FIR luminosity) and verified 
it to be robust. We investigated in detail the FIR emission of 
PACS detected and undetected sources. The average offset be- 
tween the AGN hosts and inactive galaxies is not attributable 
to larger SFRs in rare very extreme starburst events in AGN 
hosts than in inactive galaxies. It is instead due to a combination 
of modestly brighter stacked FIR emission in faint FIR sources 
(i.e. undetected by PACS) and higher PACS detection rate among 
AGN hosts compared to inactive sources, the latter feature being 
the main driver of the SF enhancement. In other words, AGNs 
are more likely to be hosted by star forming galaxies with a SFR 
above the PACS detection limit. 

The two GOODS fields, mainly characterized by low-Lx 
AGNs, and COSMOS, dominated by the brighter sources, show 
different levels of enhancement. We demonstrated that the SF 
offset is highly dependent on the X-ray luminosity: it is equal to 
~ 0.26 dex at ~ 3<x confidence level at log(Lx[erg s -1 ]) < 43.5, 
and 0.56 dex at > 10cr confidence level for more luminous AGN 
hosts. The different levels of measured FIR enhancement as a 
function of X-ray luminosity support the twofold A GN growth 
picture suggested by our previous results presented in lShao et al.l 
d2010l) . although different interpret ations cannot be r uled out 
with the present data. According to IShao et aT 1 (1201 Oh . low-Lv- 
AGN hosts undergo secular evolution, with instantaneous AGN 
feeding not strongly correlated with the total instantaneous SF 
activity in the host galaxy, and no SF offset is expected com- 
pared to inactive galaxies; luminous AGNs, on the other hand, 
co-evolve with their hosts, possibly through (major) merger in- 
teractions, and their host's SF activity is correlated with the AGN 
accretion rate. 

We suggested two hypotheses to explain the FIR enhance- 
ment observed at low X-ray luminosities. The first possibility 
ascribes the enhancement to a subset of FIR bright AGN hosts, 
and assumes that, in a major merger scenario, accretion and SF 
histories are not well synchronized, even if causally linked; a 
plausible example of this is a delayed AGN fuelling with re- 
spect to the onset of SF episodes after gas inflow, but with a long 
or uneven rise time for AGN accretion, compared to the star- 
burst timescale. A second explanation involves a larger number 
of AGN hosts, whose properties are similar to those of main se- 
quence galaxies, and it implies that a certain level of correlation 
between total instantaneous SF and accretion is also induced by 
smaller scale (non-merger) mechanisms; for example, clump mi- 
gration in massive turbulent disks drive instabilities which are re- 
sponsible for both SF and AGN accretion. Although the present 
data do not allow us to distinguish among these two scenarios 
with certainty, similarities in the SFR distributions of bright FIR 
AGN hosts and non-AGNs seem to prefer the second option. 

Future studies, especially careful morphological studies (e.g. 
in the rest-frame optical with WFC3) of luminous AGNs at these 
redshifts, will directly probe the relationship between AGN hosts 
and galaxy interactions and shed light on the process of AGN 
growth. Moreover, they will permit a discrimination to be made 



between the two scenarios suggested to explain the modest level 
of SF enhancement measured in low-L^ AGN hosts. 
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Appendix A: The main contributor to the SF 
enhancement 

As we demonstrated in Sect. 14.11 the SF enhancement in AGN 
hosts compared to inactive galaxies of similar stellar mass is due 
to a combination of two effects: modestly brighter FIR stacked 
emission in faint FIR sources (i.e. undetected by PACS) and 
higher PACS detection rate among AGN hosts. We run a simple 
calculation to understand what is the dominant effect to produce 
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the offset observed between the average vL v (60/im) of the two 
populations. 

Throughout the paper, the average vL v (60//m) in each red- 
shift and mass bin is defined as 

/ r ^ V4(«V«0 + N STACK X vL S / ACK (60fXm) 

(vL v (60/mi)) = — (A.l) 

To simplify the calculation, we assumed constant 
vL v (60yum) distributions for the detected sources, such as 
the value of each vL' v (60j-im) is equal to the average of their 
distribution. The offset between AGN hosts and non-AGNs can 
therefore be approximated by 

A(vLy(60fim)) = (A.2) 
= log(vL v (60jj.m)) AGN - log(vL v (60/im)) non _AGiv 

where 

{vLy(60fxm)) Y = (A3) 



We adopted the following assumptions which can be con- 
sidered as fairly representative of the real distributions observed 
with our dataset (Sect. l4TI ): 



- (vL° ET (60fim)) ACN = 10 45 erg s" 1 ; 

- {vL° ET {60nm))„ on -ACN = (vL° ET (60[im)) AGN : FIR bright 
AGN hosts are not enhanced in vL v (60fim) with respect to 
inactive galaxies (as suggested by the Kolmogorov Smirnov 
test); 

- vL STACK (60fxm) AG N = lO^ergs- 1 ; 

- vl) v TACK {6^m) mn - AGN = (0.4 + 1) x vL s v TACK (60fim) AGN : 
faint (i.e. undetected by PACS) AGN hosts have a modestly 
larger stacked FIR emission than inactive galaxies (see Fig. 

EJ; 

- the average detection rate (Ndet,agn /N to r,ACN) for the 
AGN population is 20% (10%), indicative of the GOODS 
(COSMOS) fields; 

- N DE T,non-ACN = (0.1 ^ 1) x Nqetagn- AGN hosts have larger 
PACS detection rates compared to the control galaxies (see 
Fig.0. 

Finally, N ST ack,y = N, otJ - N DETJ . 

We computed the values of A(vL v (60yum)) as a function 
of vL s v TACK (60[im)„ o „- ACN and N DE T,non-AGN- We first fixed 
the detection rates for the two populations to the same value 
and set vL s v TACK (60fim) mm - AC N as a free parameter: we ob- 
tained A(vL v (60jum)) < 0.08 (0.14) dex. Afterwards, we fixed 
yj^STACK (6Qjjjri) mn _ AGN to the value assumed for AGN hosts and 
allowed the detection rate for non-AGNs to vary in the allowed 
range: we found a maximum offset of 0.25 (0. 17) dex. 

Although very simple, this simulation tells that the main 
driver of the measured SF enhancement in GOODS is the higher 
PACS detection rate. Indeed, for a given detection rate, the linear 
average is highly dominated by the brightest sources, with PACS 
undetected sources only contributing in a minor fashion to the 
average offset. In the shallower COSMOS field, given the larger 
fraction of PACS undetected sources, the brighter FIR stacked 
emission and the higher PACS detection rate among AGNs can 
give comparable contributions to the average vL v (60jum) offset 
between the two populations. 
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Fig. 6. vL,,(60 / um) versus stellar mass in the different redshift bins, for the low-Lx subsample (top) and for the total one (bottom). 
The GOODS fields are represented in the upper panels and COSMOS in the lower ones. Coloured solid circles and black open 
boxes show AGN hosts and inactive galaxies, respectively. Small symbols refer to individual PACS detections, while large ones 
indicate the stacked values (for sources undetected in at least one PACS band). Error bars on the stacked values are computed 
by bootstrapping. Stellar masses for stacked sources reflect the median masses of non detected sources in each mass and redshift 
interval. In the two GOODS fields, stacked luminosities and stellar masses were averaged and weighted with the number of sources. 
Main sequence relations (grey curves) are as in Fig. [5] The inset panels show the logarithmic vL v (6Qpm) distribution for PACS 
detected AGN hosts (coloured solid histograms) and inactive galaxies (black plain ones). 
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Fig. 7. PACS detection rate in the different redshift bins, increasing from left to right, for the low-Lx subsample (left) and for the 
total one (right). The GOODS fields are represented in the upper panels and COSMOS in the lower ones. AGN hosts are showed as 
coloured solid lines and circles, inactive galaxies are represented by black dotted lines and open boxes. Icr uncertainties are derived 
through binomial statistics. Masses indicate the centre of each mass bin. Upper limits indicate cases where none of the galaxies is 
detected by PACS. Detection in all PACS bands is required. 
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Fig. 8. Same as Fig. [5] for low luminosity AGNs (see text) on the left and high luminosity AGNs on the right. The inactive galaxies 
are slightly changed, since they are randomly extracted from the total control sample to be mass-matched to the corresponding AGN 
subsample. 
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